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Abstract. The authors study a generalized notion of null geodesic defined by 
the Legendrian dynamics of a regular conical subbundle of the tangent bundle 
on a manifold. A natural extension of the Weyl tensor is shown to exist, and 
to depend only on this conical subbundle. Given a suitable defining function 
of the conical bundle, the Raychaudhuri— Sachs equations of general relativity 
continue to hold, and give rise to the same phenomenon of covergence of null 
geodesies in regions of positive energy that underlies the theory of gravitation. 

1. Introduction 

Soon after the development of the theory of gravitation by Albert Einstein, 
Hermann Weyl tried to do away with lengths and instead to base the theory entirely 
on angles (Afriat pQ, 0, [3])- In present-day language, Weyl wanted a conformally 
invariant theory, so a theory invariant under the transformations g(x) —¥ u>(x) 2 g(x), 
where g(x) is the Lorentzian metric of the theory and w(x), the conformal factor, is a 
positive function of the space-time point x. Indeed the purely gravitational degrees 
of freedom of space-time are represented by the Weyl tensor (Weyl [29]), which is 
conformally invariant (Schouten ). However Weyl's theory did not agree with 
the experiment, as was quickly pointed out by Einstein (7J. Weyl eventually recast 
and revived his theory, by turning away from gravity and instead developing gauge 
theory. In the meantime, the attempt to understand the conformal properties of 
gravity has led to much progress and has been a constant preoccupation of many 
researchers. 

The present work tries to reconcile the opposing philosophies of Einstein and 
Weyl. At the same time, we effectively generalize the theory of Einstein. The 
reconciliation is achieved by moving the entire theory into the tangent bundle of 
space-time, where the focus is on the bundle of null directions and the associated 
null geodesic foliation of that bundle, both conformally invariant constructs. In 
the Einstein theory, the bundle is defined by the vanishing of the homogeneous 
quadratic form in the velocities, g(x)(v,v) — 0, where g(x) is the metric and v 
is a tangent vector at a point x of the space-time. We generalize by allowing the 
null cone to be given by a general homogeneous function of the velocities, sub- 
ject only to the genericity condition that its Hessian with respect to the velocities 
be non-degenerate (and Lorentzian in the case of space-time). We call this struc- 
ture a causal geometry. This is somewhat similar in spirit to the construction of 
Finsler [9] geometries, although the null geodesies associated to causal geometries 
are intrinsically non-variational. 

A key point is that the allowed conformal rescalings are vastly generalized, the 
function u being replaced by a homogeneous function on the tangent bundle, so 
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a function of 2n — 1 free variables, as opposed to the standard n variables, for a 
space-time of n > 3 dimensions. So the conformal transformations are on a more 
even footing with the metric, as compared with the standard theory, where there is 
one function representing the conformal transformation and 2~ 1 n(n + 1) functions 
encoding the metric. Our first main result is that there is a natural generalization 
of the Weyl tensor, which is proven to be invariant under the enlarged class of 
conformal transformations. 

When extrapolating from an established physical theory, one wants to preserve 
as much of the structure of the old theory as possible. One motivation for going 
beyond the Einstein theory is the inevitable presence of singularities in the theory, 
as first brilliantly proved by Sir Roger Penrose [22] , [20], [21]. One might wish to 
construct a new theory free of singularities. However, at least classically, the intu- 
ition behind the Penrose theorem is compelling and depends only on the attractive 
nature of the gravitational interaction and very little on the details of the theory. 
We take this intuition as being vital to the generalized theory, so we wish to general- 
ize the Penrose singularity theorem to our case. Examination of the proof given by 
Penrose shows that apart from general causal properties described in Kronheimer 
and Penrose |17) and Geroch, Kronheimer, and Penrose [11) . which do not depend 
on the null cones being quadratic in the velocities, the only other ingredient needed 
for the proof to go through is apparently the Raychaudhuri-Sachs effect which pre- 
dicts the existence of conjugate points for congruences of null geodesies, given that 
a local positive energy condition holds and that there is a point of the congruence 
where the divergence is negative (Raychaudhuri [24], [23], Sachs [25]). Our second 
main result is that the Raychauhuri-Sachs theorem naturally generalizes to the new 
context, as does the Raychauhuri-Sachs effect in the Lorentzian case, subject to a 
natural generalization of the local positive energy condition, so this main ingredient 
of the singularity theorem goes through. 

While so far we are presenting the theory as a generalized theory of gravity, 
we see applications in many other areas of mathematics and physics, for example 
in the theory of elasticity. In our case the motivation for constructing the theory 
came from two areas studied by us: neither of these areas is concerned with gener- 
alizing the Einstein theory. Both involve the construction of a metric that is once 
degenerate and yet is not generally invariant in the degenerate direction: crudely 
speaking, a metric of the form gij (x, t)dx l dx 3 , where the (invertible) coefficient ma- 
trix gij(x,t) in general depends non-trivially on the parameter t. The question is 
what to do with this metric? 

The first example comes from the theory of third-order ordinary differential 
equations (in general non-linear) under contact equivalence |14) . We may write 
such a third-order differential equation in terms of the vanishing of an ideal of one- 
forms in four variables: {dy — pdx,dp — qdx,dq — F{x,y,p,q)dx\ where p = y', 
q = y" and the differential equation is y'" — F(x, y, y', y"). Here the prime denotes 
differentiation with respect to the variable x. 

We may pass to the three-dimensional space of solutions of the differential equa- 
tion, S. Two points of S are defined to be incident if the corresponding solutions, 
regarded as curves in the (x, j/,p)-space, meet and are mutually tangent. This in- 
cidence condition defines the null cones of an ordinary conformal structure on the 
space S, provided that a certain contact invariant of the differential equation, the 
Wiinschmann |31] invariant, W, vanishes identically. The simplest example with 
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W = 0, is the trivial equation y'" — 0, with general solution y = sx 2 + 2tx + u, 
where the null cone is that of a flat Minkowski space with conformal structure given 
by dsdu — dt 2 . 

Elie Cartan [3] and later Shiing-Shen Chern [5] studied the space, T, with co- 
ordinates (x,y,p,q). T carries a canonical direction field, V — d x + pd y + qd p + 
F(x, y,p, q)d q , such that the quotient of T by V is the space S. Chern showed that 
T carries a once-degenerate conformal metric, which is killed by V , such that it is 
also invariant under V (up to scale), so passes down to S, if and only if W = 0. 
In the case that W ^ 0, the present authors showed that one could use the Chern 
metric to construct on the space S a null cone structure by the method of envelopes 
and which reduces to the standard null cone structure in the case W = 0. This null 
cone structure is exactly that given by the incidence condition. Thus the causal 
geometry is natural for this case and one wants to develop an analogue of the 
usual connection theory which applies in this case. Our theory does this, although, 
ironically, the Weyl curvature vanishes identically, as it does in the standard case 
of W = 0, for dimensional reasons. 

The simplest example with W ^ is the differential equation y'" — y". Its 
solutions are y = se x + tx + u, where the parameters (s, t, u) are global co-ordinates 
for S. The incidence conditions arc = dy — pdx = e x ds + xdt + du and = 
dp — qdx = e x ds + dt. Eliminating x between these equations gives the causal 
null cone in the form e 1- "^ + ^| = 0. This is well defined and has non-singular 
hessian with respect to the variables (ds,dt,du), provided only that dt ^ 0. It is 
dramatically more complicated than the case of W = 0. 

A remaining issue is to relate this work to that of the second author and Pawel 
Nurowski, who built a canonical conformal structure in six dimensions that en- 
codes the geometry of the third-order equation and which, when W = reduces 
to a conformal structure of the type first given by Charles Feffermann [8]. When 
W 7^ 0, one needs a generalized Fefferman structure, applied to general parabolic 
geometries, as shown by Hammerl and Sagerschnig |13j . 

The second example comes from the twistor theory of four-dimensional real 
curved space-time. One considers the co-spin bundle of the space-time. Points 
of the bundle are written (x, ~ka')i where tta> is a two-component dual primed 
spinor and we use the abstract spinor and tensor indices of Penrose. The co-spin 
bundle is a C 2 bundle over space-time, so has total real dimension eight. Briefly 
the co-spin bundle carries an Ehresmann connection, written diVA 1 , with complex 
conjugate (Ma, which encodes the Levi-Civita spin connection of the space-time. 
Then the co-spin bundle carries a canonical degenerate real metric, J 7 , called the 
Fefferman metric [27] : 

J- = i6° '(w ' Ad'n a 1 ~ ^a'cMa)- 

Here 9 a is the canonical one-form of the space-time. It is not difficult to show 
that T depends only on the conformal structure of the space-time. It is actually 
degenerate in two directions. The first is that of the vector field generating the 
scaling transformations tta' — > tit A' > where t is real and positive. Then T is invariant 
up to scale under these transformations, so one may factor out by identifying tta' ^ 
with t-KA' for any positive t, giving an § 3 bundle over space-time on which T 
gives a degenerate conformal structure. The remaining degeneracy is that of the 
null geodesic spray, V = tt a Tf A d a: where d a is the system of horizontal vector 
fields of the Ehresmann connection on the spin bundle, representing the Levi-Civita 
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connection. Then the conformal structure J- is invariant under the null geodesic 
spray, if and only if the Weyl curvature vanishes, if and only if the space-time is 
conformally flat and in that case J- passes down to the quotient space as a non- 
degenerate conformal structure, which is actually a standard Fefferman conformal 
structure for the space. So the vanishing of the Weyl curvature gives the analogue 
of the condition of vanishing Wimschmann condition in the case of third-order 
differential equations. In the case that the Weyl curvature does not vanish, J- 
rcstricts to give the Fefferman conformal structure of a canonical C1Z structure for 
the projective spin bundle over any hypersurface in the space-time, which is its 
defining property. As in the case of third order ordinary differential equations, it is 
not clear what to do with T in the general case. One strategy is to use the envelope 
approach as in the third-order case. Then one obtains a causal geometry on the 
space of null geodesies. 

There is only one non-conformally flat case where the relevant calculations have 
been carried out in detail, that of the Kapadia [15] metric: 

g = dudv — dx 2 — u~ dy . 

Here the co-ordinates (u,v,x,y) E M 4 with u > 0. The metric in null and is con- 
formal to vacuum. The null geodesies and the null cones are explicitly computable, 
with the result that the null cone of a point (uq, vq,xq, yo) is given by the formula: 

n / v ^ ( ^2 2 (y-Vo) 2 

0= {u-Uo){V-V ) -[X-Xq) ; . 

u + u 

After some calculation, one can put the Fefferman conformal structure in the form: 

T = 2dO{dV + 2PdX + AQdY) + u^^dX(dY + 3u 2 dQ) - u~^dP{3dY + u 2 dQ). 

In these co-ordinates, the degeneracy direction is that of the vector field V = J^. 
Differentiating with respect to u, we get: 

2 5 

--ui^yT = (dX - udp)(dY - u 2 dQ). 
o 

The causal structure for the quotient by V is obtained by eliminating the variable 
u between the relations J- = I£yT = 0, so splits into two parts: 

• For the first, with dX — udP, we have, after a rescaling, the causal structure 
M: 

M = dd(dV + 2PdX + 4QdY) (—) - dYdP + dXdQ—. 

\ dP J dP 

Then A/"i is well defined and non-degenerate, provided > 0. 

dP 

• For the second, with dY = u 2 dQ, we have, after a rescaling, the causal 



structure Af' 



2- 



' dY\ 3 I dY 

A/2 = 2d6(dV + 2PdX + AQdY) ( — + AdYdX - AdYdP 



dQ J \dQ 
dY 

Then A/b is well defined and non-degenerate, provided — — > 0. 

dQ 

In both cases the null geodesic equations are completely integrable by quadratures. 
In both cases the tidal curvature may be calculated explicitly; in each case, when 
written out, the curvature involves some 533 terms. 
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2. Preliminaries 

2.1. Secondary bundle structure on the tangent bundle. Let M be a smooth 
manifold of dimension n > 2. The tangent bundle TM of M consists of pairs (x, v) 
with x £ M and w £ TM Xl the tangent space to M at x. The bundle projection 
ittm ■ TM ^ M is defined by ^tm(x,v) — x. The double tangent bundle TTM 
is the tangent bundle of the tangent bundle, and consists of triples (x, v, w) where 
(x,v) £ TM and w £ TTM {x ^ v) . The bundle projection tt T tm ■ TTM -> TM is 
defined by 7Tttm(x, v, u>) = (x, v) £ TM. 

In local coordinates x % of M, there are induced linear coordinates v % in each fiber 
of TM, defined by 

i/ ^ d 
V = v{v) dx^ 

Then TTM also carries fiber coordinates in the 2n-dimensional space TTM( X;U ), 
denoted by £*, i/*, defined at w £ TTM( X; „) by 



ax* 



7/ N ^ 

[x,v) 0V 



(x,v) 



Apart from the bundle projection ttttm '■ TTM — > TM on the second tangent 
bundle, there is also another natural projection given by the differential (Tktm ■ 
TTM — > TM. In the local coordinates described above, 

d d 

dTTTMTT- = -K— 

ox ox 1 
d 

dlTTMTT- = 0. 

The kernel of (Itttm is called the vertical subbundle, and is denoted by VTM. 
There is a natural isomorphism between VTM and the pullback bundle tt^^TM, 
given as follows. Let x 6 M and v,w £ TM X . The one-parameter group L w (s) : 
(x, v) n> (x, v + sw) as s £ M varies, is a well-defined one-parameter group of 
diffcomorphisms of TM X to itself. Denote the generator of this one-parameter 
group by \ x , v )(w) = L' w (0). Then A^^) : TM X ->■ VTM^^y This is a linear 
isomorphism for each fixed (x, v) £ TM, and it depends smoothly on (x, v). So it is 
an isomorphism A : n^^TM — > VTM of vector bundles over TM. In coordinates, 

A— — 

dx % dv % 

Since tx^\/TM = TTM '/VTM ', it is convenient to compose A with the quotient 
map q : TTM -> TTM /VTM to obtain A = A o q : TTM -)• FTM. Then the 
image and kernel of A are both the vertical bundle VTM. As a tensor, A can be 
identified with a section of V a TM ® VTM, In coordinates, 

A— — A— 

dx % dv l ' 9u l 

and, as a tensor, A = dx l ® J^. 

Let A be a vector in TTM. Define a differential operator D x : C°°(TM) -> 
C°°(TM) by 

D x {f)=J? xx f. 
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In local coordinates, if X = f gfr+^^r, then D x (f) = C J£- Let D : C°°{TM) -> 
rrM(? 1 *rAf) be the one-form valued operator 

(Df){X) = D x f. 

If X, y G IVm (7r,p^ TM) are * wo vec tor fields that are lifts of vector fields from 
M, then 

D\y = D X D Y = D Y D X , 
and D x Y (f) depends bilinearly on X, Y. Commutativity follows from the commu- 
tativity of the one parameter groups L x and Ly defined previously. 

Let TM' be the tangent bundle with the zero section removed, and tttm' = 
ktm\tm> ■ TM' -> M the induced projection of TM' onto M. Let SM be the 
space of oriented one-dimensional linear subspaces of TM. Let S s : TM' — > TM' 
be the scaling 5 s (x, v) — (x, sv) for s > 0, and let H — 4-6 s \ s —x be the homogeneity 
vector field. This defines a group action of (0, oo) on TM', and SM is the quotient 
bundle of TM' by the group. Let itsm '■ SM — > M be the projection onto M. 
There is a factorization nrw = ^SM ° c where <r : TM' — >• SM is the quotient 
mapping. 

2.2. Frolicher Nijenhuis bracket. Let X be a smooth manifold and f2(X) = 
© r f2 r (X) be the graded algebra of smooth differential forms on X. A derivation of 
degree k of f2(X) is a real linear map £) : ft(X) — > fl(X) such that 

• D : fl r (X) -> ft' r+fe (X) 

• For any a G O a (X) and /3 G n b (X), D(a A /3) = (Da) A /3 + (-l) fca a A D(3 
Let Derfc(r2(X)) be the space of derivations of degree k of Q(X), and let Der(f2(X)) = 

©fcez Derfc(fi(X)) be the graded vector space of all derivations; this supports the 
structure of a graded Lie algebra, where the bracket of homogeneous elements 
K G Der k (n(X)),L G Dcr e (n(X)) is defined by 

[K,L] =KoL-(-l) ke LoK. 

Extending by bilinearity to all of Der(f2(X)), the resulting bracket is easily seen to 
define a graded Lie algebra: 

• The bracket is graded anticommutative: 

[K,L] = -(-l) ke [L,K] 

for K G Der k (Q(X)),L G Dcr e (n(X)) 

• The bracket satisfies the graded Jacobi identity: 

(-l)'V, [K,L]} + (-l) k i[K, [L,J]} + (-l) ek [L, [J,K]} = 0. 

If v € T X {TX) is a vector field, then the insertion operator i v : ft(X) —> ft(X) 
is a derivation of degree —1. The insertion operator extends to an operator i^ G 
Derfc_i(f2(X)), by defining i^t&v = ui v for K G fl k (X, TX) extending by linearity. 

Definition 1. Let K G £l k (X, TX). Define the Lie derivative along K by 
%K = [ik, d] =i K oal+ (-l) k d o i K . 
The following is proven in |16j : 
Theorem 1. Any derivation D G Derfc(ri) can be decomposed uniquely as 

D = Jz?a' + i L 
for some K G fl k (X : TX) and L G fl k+1 (X,TX). 
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Definition 2. For Ktl k (X ,T X) and L e Vt l (X,TX) define the Nijenhuis-Richardson 
bracket [K, L] A by 

Define the Frolicher-Nijenhuis bracket [if, L] &y 

The algebraic properties of the curvature and related quantities are most easily 
expressed using the Nijenhuis-Richardson bracket [18] and the Frolicher-Nijenhuis 
bracket [TO] . 

Lemma 1. [JS?jc,ii,] = - (-l) M i^ L K 

Proo/. □ 
Lemma 2. 

[K, [L U L 2 ] A ] - [[if, LJ, i 2 ] A + (-1)** [Lx, [K, L 2 ]] A - 

- ((-l)**[i Ll tf,L 2 ] - (-l^+W^K,^)) 
Proof. D 

The following is also useful: 
Lemma 3. IfK,L E Q}{X,TX), then 

[K,L](A,B) = [KA, LB] - [KB, LA] - L([KA, B] - [KB, A]) 

- K([LA, B] - [LB, A]) + (KL + LK)[A, B). 

3. Null geodesic dynamics 

3.1. Legendrian dynamics. The dynamics is specified by a smooth hypersurface 
<£ (of dimension 2n — 2) in SM, with 7f(^) = M, which has the property that for 
any x € M, the intersection Sf R SM X is a smooth submanifold of the sphere SM X 
of dimension n — 2. Note that the space Jf? = o~ lc S ', a smooth hypersurface in 
TM' invariant under the scaling (x,v) —> (x,tv), also specifies the dynamics. 

If (x,v) G Jtif, the vertical tangent space VM[ x ,v) to Jf? at (x,v) is the in- 
tersection of the tangent space to Jff at (x,v) with the vertical space VTM^ xv y 

So VM^ x ,v) has dimension n — 1. The image of VJ^ X V ^ under the map A^ K \ : 
7TM(j. „\ — > TM X is then a subspace of TM X , also of dimension n — 1. There 
is then a unique maximal subspace of TJf?( x ^ that projects down under the map 

d/KTM' to A^ v ^(VJ^f( XtV -j). This subspace is a codimension one distribution within 
TJff, and therefore defines a distribution of hyperplanes on Jf?; 



This entire construction is invariant under the scalar homothety S, and so Ajp 
descends to a distribution of Ag? on §f as well. 

Definition 3. ^4 contact symmetry of a distribution A on a manifold X is a one- 
parameter local group of diffeomorphisms of X that preserves A and whose genera- 
tors are everywhere tangent to A. 

Definition 4. The dynamics o/Sf is the space of contact symmetries of Acg . 
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3.2. Lagrangian approach to the dynamics. Let <#C have local defining equa- 
tion G(x, v) = where G is a smooth function defined over an open set U of TM' 
that is invariant under 5 satisfying: 

• DG ^ throughout U 

• G homogeneous of some real degree k: G(x,tv) = t k G(x,v) for all t > 
and all (x, v) G U. For convenience, we shall henceforth assume that k ^ 1. 

There is a bilinear form g/j on -k^ m ,TM defined for vector fields X and Y that 
lift vector fields on M by 

g h (X,Y) = D% Y G. 

The definition is independent of the choice of lift of X and Y, and it is bihomoge- 
neous under rescalings X — > (Tr^ M ,a)X and Y — > (iT^ M ,b)Y where a, b are functions 
on M, and so it gives rise to a bilinear form. The h here stands for "horizontal" , a 
reflection of the fact that gh is a section of V°TM ® V°TM. Applying A yields 
a bilinear form g v in V*TM <g> V'TM: 

5 „(x,y) = .g ft (A" 1 (x),A" 1 (r)). 

Here the subscript w means "vertical" , since g v is a bilinear form on VTM. In 
coordinates, 

<9 2 G 

5d = dvl Q v j dyt ® = 9tjdv l O dw 3 

flh = g~[Q~j dxt ® dx3 = <8> rfa; 3 ■ 

Lemma 4. Lef a e T TM ,(T*TM') be the differential form a = DG. Then, on 
restricting to ,yiC , the distribution Aj? is the annihilator of a in TM'. 

Proof. By the assumption that DG ^ 0, the image of TJf under a is always 
one-dimensional, and so the annihilator of a is a distribution of hyperplanes in 
TJf. Suppose a(X) = for X e TJf. Then, by definition of the D operator, 
XX j dG = 0. So XX e VM. That is, X e X~ 1 {VM') as required. □ 

Lemma 5. g h (X,Y) = 2da(XX,Y) 

Proof. Both sides vanish if either X or Y is vertical, so it is sufficient to establish 
the lemma under the additional assumption that X and Y are lifts of vector fields 
from M. Since a(XX) = 0, 

2da(XX,Y) = (& xx a)(Y) = XX(a(Y)) - a([XX,Y}) 

= D 2 X Y G - a([XX, Y}) - g h (X, Y) - a([XX, Y}) 

But if X and Y are lifts of vector fields, then [XX, Y] is vertical, and so a([XX, Y]) = 
0. □ 

Assume henceforth that the bilinear form g v is nondegenerate. This assumption 
is justified in part by 

Lemma 6. The bilinear form g v is nondegenerate if and only if da is a symplectic 
form on a neighborhood of in TM'. 

Proof. The subspace VTM' is an isotropic space for da. Choose a complementary 
space HTM' in TTM'. Then da induces a bilinear form on VTM' x HTM' and 
da(X,Y) = 2g v (X,XY). □ 
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In coordinates, 



a = ——dx 



Set pi — dG/dv 1 . By the nondegeneracy of g, the Jacobian matrix dpi/dv^ is 
nonsingular, and so this defines a new set of (local) coordinates on TM'. In the 
new coordinates, 

a = pidx 1 . 

These are the "canonical coordinates" for the dynamical system. 

The symplectic form da allows us to define the Poisson bracket of two functions 
fi , fi ( m a neighborhood of J^) by 

{fi,f2} = (da)- 1 (df 1 ,df 2 ). 

This satisfies the usual rules: 

• {/l,M = -{/2,/l} 

• {/i, 1/2, / 3 }} + {/ 2 , {/s, /i}} + {/a, {/i, / 2 }} = 

• c l = if c is constant 

• {/l,/2 + /3} = {/l,M + {/l,/3} 

• {/l,/2/3} = {/l,/2}/3 + {/l,/ 3 }/2 

The last three properties imply that the operator } : / 2 >->• {/i,/ 2 } is a 

derivation on smooth functions, and therefore corresponds to a vector field on M. 
In the canonical coordinates, 

dpi dx l dx l dpi 

As we are interested in the intrinsic geometry of Jif, we shall consider the pull- 
back of a to Jff. 

Lemma 7. When pulled back to Jff, a has Darboux rank 2n — 3: 

a A {da) n - 2 ^ 0, (da) 11 - 1 = 0. 

Proof. The fibers of J^f — s> M are n — 1 dimensional, and the bilinear form (/„ 
on VJif is annihilated by the generators of scaling up the fiber. So on VJif, g v 
has rank n — 2. By the argument in the previous lemma, {da) 11 ^ 1 — 0. However, 
applying the previous argument to Aj^> = a , and choosing a complement for this 
in TSf gives a A (da) n ~ 2 ^0. □ 

The first main result uses the Darboux theorem: 

Lemma 8. Let M be a manifold of dimension 2n—l and a a one- form of Darboux 
rank 2r — 1 . Then the space of vector fields X such that 

(1) Xja = 0, aAJz?xa = 

forms an integrable distribution of rank 2(n — r). 

Proof. Using the usual version of Darboux' theorem, there exists a coordinate sys- 
tem x, y 1 , . . . , j/ n_1 ,pi, • • • ,Pn-i on M such that 



a = dx + y^Pi oh/, 
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with the last (n — r) p's and y's not participitating. Hence the vector fields d/dpi 
and d/dy l for i = r, ...,n — 1 form an intcgrable distribution of rank 2(n — r) 
satisfying ([1]). 

Now, note that any X satisfying (p} must also satisfy 

X^aAidaf- 1 ) = 0. 

But a A {daY^ 1 = ±dx A dpi A • • • A dp r -i A dy 1 A • • • A c?y r ~ 1 is annihilated by X if 
and only if X is a linear combination of d/dpt and d/dy 1 for i = r, . . . ,n — 1. □ 

Theorem 2. TTie dynamical vector fields on J^f are spanned as a C°° module by H 
(the generator of the scaling symmetry of M') and the vector field V = (k — l){G, — } 
restricted to 

The particular normalization of V ensures that it defines a spray; see Lemma [S] 
below. 

Proof. By Lemma there are exactly two linearly independent dynamical vector 
fields at every point. Note that V is tangent to since V{G) — (k — 1){G, G} = 0, 
and H is tangent to Jif since Jif is invariant under the scaling action. These are 
linearly independent, since V(ir*%,f) — (k — l){G, 7r*^./} is nonzero for some smooth 
function / on M, but H{-K*^f) = for all such /. 

Now, note that H satisfies H _ia — (since XH = 0). If X is a lift of a vector field 
on M, then 6 S XX = s _1 AX. Differentiating gives ££h(XX) = —XX. For such a 
vector field X, [H,X] is vertical and so X[H,X] = 0. It follows that a([H,X\) = 0, 
and therefore 

(& H a){X) = ^ H {XXjdG) -a([H,X}) 

= -XX j dG + XX j ^ H dG 

= {k- l)XXjdG =(k- l)a(X). 

Finally, in a neighborhood of Jtf, V is characterized by 

Vjda = -(k-T)dG. 

Pulling back to G — gives V-ida = 0. By the previous calculation, H j da = 
(k — l)a. Hence 

7j a = (k - H_,da = -(k- l^Hj Vj da 

= H j dG = kG 

which also vanishes on Jff. □ 
In coordinates, the dynamical vector fields are 

^ V dv l 

V -v l —+u l — u l - dG v* 

dx l dv l ' y dxi dx l dvi 

A direct derivation of the previous results in coordinates is given in an appendix. 

The integral curve of the vector field V through a point (x, v) projects to a curve 

in M whose initial velocity is v. That is, V is a semispray. 

Lemma 9. V is a spray: 
• [H, V] = V 



H = XV 
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Proof. The calculations in the proof of the preceding lemma give 
[H, V]j da = ^ h {Va da) - Vj(^ H da) 

= k(Vj da) - (jfe - 1)Vj da = Va da, 

so [H, V] = V. 

For the second property, the definition of gh implies 

da(X,XY) = ^g h (X,Y) = da(Y,\X) = -da{XX,Y). 

In particular, with X = V, 

da(XV, Y) = -da(V, XY) = (k - l)dG(XY) = (k - l)a(Y) = da(H, Y). 
This is true for all Y and so H = XV by nondegeneracy of da. □ 

Passing down to the sphere bundle SM, only the dynamical vector field V sur- 
vives, up to an overall positive scale, since H is in the kernel of da : TTM — > TSM . 
This gives a foliation of Sf by the dynamical curves, the (maximally extended) tra- 
jectories of V. These dynamical curves are called null geodesies. The space of null 
geodesies, denoted by Jf , has dimension 2n — 3. The distribution keg is Lie derived 
along the dynamical vector fields, and so descends to a codimension one distribution 
on JV . This distribution is a contact structure since the relation a& f\(day ) n ~ 2 ^ 0, 
valid for any nonzero a<g in the annihilator of Ay, also descends to the quotient. 

The null geodesies are naturally oriented, since at each point p of <S ', the vector 
field V descends to a ray through the origin in Tj0 ', which is oriented. The bundle 
'S is time oriented if and only if the space of oriented null geodesies is the disjoint 
union of two components, ,jV — J/ + U , such that the oriented null geodesic 
through (x, v) lies in ,jV ± if and only if the oriented null geodesic through (x, —v) 
lies in ,jV^ . Then the elements of are called future oriented and the elements 
of ,jV~ arc called past oriented. 

On ffl , the integral curves of V are called affinely parametrized null geodesies. 
These carry a natural parametrization up to a translation, since they are the integral 
curves of a single vector field. This natural parametrization requires having a 
particular defining function G for Jff, although the definition of (unparametrized, 
oriented) null geodesies on £f does not. 

4. Ehresmann connection 

4.1. The Ehresmann connection on Jf. The purpose of this section is to es- 
tablish the following: 

Theorem 3. There exists a unique operator P : TTM' — > VTM' satisfying for all 
X, Y e TTM: 

(1) (J? v g h )(X, Y) = g v (PX, XY) + g v (PY, XX) 

(2) 2da(X, Y) = g v (PX, XY) - g v (PY, XX) 

This operator defines an Ehresmann connection on TM' , meaning that it has max- 
imal rank and satisfies P = P o P. Furthermore, 

(2) P=\ (Id TTM ' +J^VA) . 

The restriction of P to is also an Ehresmann connection on : P(TJf?) = 
VJf. 
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(Recall that A is a section of V°TM'®VTM' , g h is a section of V°TM'®V°TM', 
and g v is a section of V*TM' ® V*TM'.) 

The proof is broken down into several lemmas. 

Lemma 10. The Frolicher-Nijenhuis bracket of X with itself is zero: [A, A] = 0. 
Thus if X,Y eT TM >(TTM'), then 

\{[\X, Y] + [X, XY]) = [XX, XY}. 

Moreover [V, XY] = -Y (mod VTM') for all vector fields Y. 

Proof. If X,Y E T T M'(TTM'), then 

i[A, X](X,Y) = [XX, XY] - X([XX, Y] + [X, XY]) 

The right-hand side vanishes if X or Y is a section of VTM', since VTM' is an 
integrablc distribution on which A vanishes. Thus it suffices to prove that it vanishes 
if X and Y are both lifts of vector fields from M . In that case, if / is a function 
on M , then 

[XX,Y]n* TM ,f = (XX)Yir* TM ,f = (XX)w* TM ,((dir TM/ Y)f) = 0. 

Hence [XX, Y] G T TM >{VTM'); likewise [X, XY] e T TM ,{VTM'). So X([XX, Y] + 
[X, XY]) = 0. Finally, [XX, XY] = as well for X,Y lifts of vector fields on M, 
since the one-parameter groups Lx and Ly (defined in i j2.ip commute in that case. 

It remains only to show that [V, XY] = — Y (mod VTM') for all vector fields Y. 
Taking X = V in the first part gives 

X([XV,Y] + [V,XY]) = [XV,XY]. 

But Ay = H, so rearranging gives 

a[v,ay] = (jsf ff A)(y) = -y 

as claimed. □ 

Lemma 11. ^4ny operator P : TTM' — > VTM' satisfying property (1) has maximal 
rank and satisfies P o P = P. 

Proof. In view of the fact that imPc VTM' by assumption, it is enough to show 
that P(X(X)) = X(X) for all X. This then proves that the range of P is equal to 
the vertical tangent space, and that P acts as the identity on its range. Therefore 
P o P = P, and P has maximal rank. 
To prove the claimed identity, (I) gives 

{^v9h){XX, Y) = g v (PXX, XY) + g v (PY, XXX) 

= g v (PXX, XY) 

since XXX = 0. Expanding the left-hand side, 

(J? v9h )(XX, Y) - V(g h (XX, Y)) - g h ([V, XX], Y) - g h (XX, [V, Y]) 

= -g h ([V,XX],Y) 

since every vertical direction lies in the kernel of gh- But, by Lemma [HI [V, AX] = 
-X (mod VTM), and therefore 

-g h ([V,XX},Y)=g h (X,Y). 



CAUSAL GEOMETRIES, NULL GEODESICS, AND GRAVITY 



13 



Putting these together, 
(3) g h (X,Y)=g v (PXX,XY). 

This is true for all X, Y, and so PXX = XX, as claimed. □ 

Lemma 12. There exists a unique P G T T M'(T*TM' ® VTM') satisfying condi- 
tions (1) and (2). 

Proof. If (1) and (2) hold, then 

g v (PX, XY) = l - ((J? v g h )(X, Y) + da(X, Y)) . 

By the non-degeneracy of g v and the fact that A : TTM' — > VTM' has maximal 
rank, this admits at most a unique solution P(X) valid for all Y. To prove existence, 
it is enough to show that the kernel of 

Y^(J? v9h )(X,Y) + da(X,Y) 
contains the kernel of A, which is also the image of A, VTM'. So consider 
(J? v g h )(X, XY) + da(X, XY) = -g h (X, [V, XY}) + da(X, XY) 

= g h (X,Y)-g h (X,Y)=0 

where we have used the fact that V is a spray in simplifying the first term, and the 
definition of g^ in simplifying the second term. □ 

Lemma 13. The unique connection P satisfying (1) and (2) is given explicitly by 

P=\ (ld TTM , +J? V X) . 

Proof. Let P be the connection characterized by (1) and (2) and let P\~\ (Wttm' +-2?y 
We will show that Pi acts as the identity on VTM', and that ker Pi = ker P. The 
first claim follows at once from 

(^f v X)(XX) = [V, X 2 X] - X[V, XX] = XX. 

For the second claim, X G kerPi if and only if («SfyA)(X) = —X, or, equivalently, 

-X = [V,XX] - X[V,X\. 

Now X G ker P if and only if 

(J?vgh)(X,Y) = -2da(X,Y) 

for all Y G TTM'. Note 

(J? v g h )(X,Y) = V(g h (X,Yj) - g h ([V,X],Y) - g h (X, [V,Y]) 

= V(da(XX, Y)) - 2da{X[V, X],Y) — 2da{XX, [V, Y]) 

= 2da([V,XX],Y) - 2da(X[V, X],Y) 

= 2da((J? v X)(X),Y). 

If X G ker Pi, then this last display reduces to — 2da(X, Y), and sole ker P. Con- 
versely, if X G kerP, then the same calculation shows that 2da((^vX)(X), Y) = 
—2da(X,Y) for all Y, and hence (^fvX)(X) = —X by the nondegeneracy of da, 
and so X G ker Pi. □ 
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In coordinates, 

i , ; tt, i i\ d rri 1 du l , dG , d 2 G 



J ' <9u* J 2 9t>i ' ay 9xJ 0a*9uJ ' 

The horizontal lift of the coordinate vector fields d/dx 1 are 

h{d/dj) = (I - P)(9/9^) = Jjg-uiJL. 
Lemma 14. TTie operator P : TTM' -> satisfying (1) and (2) is such that 

on jr, p(tjt) = yjr. 

Proof. Let X be a vector field in that Lie commutes with V. Since AV^ = H , 
g v (P(X),H) = g v (P(X),XV) = l - (da(X, V) + (Sf Y g h ){X, V)) . 

It is sufficient prove that both terms of the right-hand side are zero. By the calcu- 
lations in the proof of Theorem [2] 

da(X, V) = (k — l)dG(X) = 

since X is tangent to Jff. Also, since X, V commute by hypothesis, 

(Sf V 9h)(X, V) = V(g h (X,V)). 

Now g h (X, V) = -2da(X, XV) = -2da(X, H) = (k - l)a(X), and thus 

V (g h {X, V)) = {k- l)V(a(X)) = (k - l)(J? v a)(X) = (k - l)dG(X) = 0. 

□ 

5. Curvature 

5.1. Tidal force. 

Theorem 4. Let P be the connection of Lemmas \ll\ and \12\ Then there exists S G 
T TM >(V°TM' ® V°TM') such that \{^g h ){X,Y) = g v (P(X),P(Y)) + S{X,Y) 
for all X, Y G TTM' . Conversely, P is the unique operator such that 

(1) (3?v9h)(X, Y) = g v (P{X),Y) + g v {X, P(Y)) for all X.Y e TTM'. 

(2) There exists S G V°TM'®V°TM' such that \{^g h ){X, Y) = g v {P(X), P(Y)) + 
S(X, Y) for all X, Y G TTM'. 

The symmetric tensor S is called the tidal force tensor. 

Proof. For the first claim, it is enough to show: 

(1) {S%g h )(\X, XY) = 2g v (XX, XY) 

(2) (3e$g h )(\X, Y) = for all X and Y in the kernel of P. 

Indeed, assuming these are both true, decomposing two vectors X = Xh + X v and 
Y = Yh + Y v into ker P and im P components, 

(&$g h )(X h +X v ,Y h + Y v ) - g v {X v ,Y v ) = (£>*g h ){X h , Y h ) 

which defines S(X, Y). 

For (1), since (& v g h ){\X, XY) = and g h (XX, Z) = g h {Z, XY) = for all Z, 

(J?*g h )(XX, XY) = -(^ v g h )([V,XX],XY) - (J? v g h )(XX, [V, XY]) 

= 2g h ([V,XX},[V,XY}) = 2g h (X,Y) = 2g v {XX, XY). 
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For (2), Y is in the kernel of P if and only if 

{Jf v g h )(Y,Z) + 2da(Y,Z) = 

for all Z. Hence 

(J? v g h )(XX,Y) = V((J?v9h)(\X,Y)) - (J? v g h )([V,XX],Y) - (J? v g h )(XX, [V,Y]) 
= V(g v (P(XX), XY)) + da(Y, [V, XX}) - g v (P(XX), X[V, Y}) 
= V(g h {X, Y)) + da(Y, [V, XX}) - g h (X, [V, Y}) 
= (J?v9h)(X, Y) + g h {[V, X], Y) + da(Y, [V, XX}) 
= g h ([V,X},Y) + 2da(Y,[V,XX}) 

since da(X, Y) = for X, Y 6 kerP. Now, 

[V, XX} = (& V \){X) + X[V,X} = (2P - ld)X + X[V,X] 

= -X + X[V,X}. 

So, continuing the above calculation gives 

(J?$g h )(XX,Y)=g h ([V,X],Y)+da(Y,X[V,X})=0 

where we have used again the fact that da(X, Y) = along with Lemma [5] 

For the converse statement, let £P be the affine space consisting of all operators 
P : TTM' ->■ VTM' satisfying 

(J?v9h)(X, Y) = g v (P(X),X(Y)) + g v (X(X), P(Y)). 

By Lemma [TTJ any such P satisfies P o P = P 1 and so defines a projection onto 
VTM' . Any such operator is completely determined by its kernel. But for the P 
satisfying Lemma [T^l it follows from the first part of the lemma that 

kerP= p| ker[(^g h )(X,-)}. 

XeVTM' 

□ 

5.2. Curvature of the connection. The curvature of the Ehresmann connection 
P is the section of V°TM' ® V°TM' ® VTM' defined by 

R(X,Y) = P[(ld-P)(X), (Id-P)(Y)]. 

Since the vertical bundle is integrable, this evaluates to 

R(X,Y) = [PX.PY] - P{[PX,Y] + [X,PY})+P[X,Y\. 

Equivalently, this can be re-expressed in terms of the Frohlicher-Nijenhuis bracket 
[IS], by 

Lemma 15. The Bianchi identity holds [P, R} = 0. That is, 

[PX, R(Y, Z)]+P[R(X, Y), Z]+R([X, Y], Z)-R([PX, Y],Z)+R([PY, X], Z)+cycKc = 0. 

Define S x e T*TM' ® VTM' to be the unique tensor such that 

g v (S x X, XY) = S(X,Y) 

for all X, Y, where 5 is the tidal force tensor of Theorem 01 The curvature deter- 
mines the tidal force, and vice versa: 
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Theorem 5. The curvature and tidal force are related by 

S X X = -R(V, X) 

for all X, Y £ TTM'. Moreover, 

[X,S X ] = ~R. 
The following lemma is of interest in its own right: 
Lemma 16. The tidal force tensor satisifies 

S(X,Y) = ~g v (P(J?$X)X, XY). 

So S x = iP(_Sf£A). Moreover, if X,Y e kerP, then 

(4) S(X,Y) = -g v (P[V,X],XY), 
Proof. Note first the operator identities 

(5) P(^VA) = P, («5fyA)A = A. 
Thus 

P(^X)XX = P[V, [5? V \)\X] - P(& V X)[V, XX] = 0. 

So it is sufficient to establish the first statement of the lemma under the addi- 
tional hypothesis that X, Y £ kcrP. In that case 

X = -{5? V X)X 

so 

P{^X)X = P[V, (JSVA)X] - P(^ V X)[V, X] 

= -P[V,X]-P[V,X] = -2P[V,X]. 

So to prove the first part of the lemma, it is enough to show Q . 

Since X, Y £ kerP, it follows by Theorem H that 2S(X,Y) = (^g h )(X,Y). 
Now, 

2S(X, Y) = (J?$g h )(X, Y) = V (J?vg h (X, Y)) - Sf v9h ([V, X],Y) — J? v g h (X, [V, Y}) 

= -2da([V,X],Y) +2da(X, [V,Y]) 
= 2V{da{X,Y)) -4da([V,X],Y) 
= -4da([V,X],Y) 

since da(X, Y) — for X, Y £ kerP by Theorem [3] Thus, applying Theorem [3] 
twice more, 

2S(X,Y) = -4da([V,X],Y) = -2(J? v9h )([V,X},Y) = - 2g v (P[V,X], XY) . 

□ 

Proof of Theorem^ For the first identity, it is enough to show that (Jz?ygh)(X, Y) — 
-2g v (R(V,X),XY) for all X, Y £ kerP. Since V £ kerP as well, R(V, X) = 
P[V,X], and so the last part of Lemma [TBI gives 

S(X,Y) = -g v (R(V,X),XY), 

as claimed. 



CAUSAL GEOMETRIES, NULL GEODESICS, AND GRAVITY 17 

Now, by the first part of the theorem, S x = —\iyR- So, by Lemma [21 
2[A,S A ] = -[A, [V,R] A ] 

= - ([[A, V], P] A + [V, [A, P]] A - [i v X, R] + [i R \, V]) . 
Now ijiX — 0, iyX = Ay = H, and 

[H,R] = [[H,P],P} + [P,[H,P]]=0 
[[A, V],P] A = -L2VA,P] A = -[2P-Id,P] A 

= -2[P, P] A + [Id, R] A = 2R + R = 3P. 

So 

(6) [A,S X ] = -(3P+[V,[A,P]] A ). 

Now, we claim that [X,R] = 0. By the graded Jacobi identity, 

[A,P] = ~[A, [P,P]] = -[P,[A,P]]. 

From P = i(Id+«5f y A), 

[A,P] = i([A,Id] + [A, [V,A]]) = i[A > [V, A]] 

since [A, Id] = by Lemma [31 The graded Jacobi identity applied once more gives 

-[A, [V,X]] + [V,[X,X]] + [X,[X,V}] = Q. 

But [A, A] = by Lemma M and [V, A] = - [A, V}. So [A, [V, A]] = 0, and therefore 
[A, R] = 0, as claimed. 

So becomes [A, S x ] = — |P as required. 

□ 

Since R is skew-symmetric in its arguments, the first part of the theorem implies 
immediately 

Corollary 1. For any X G TTM', S(V,X) = 0. 

6. CONFORMAL TRANSFORMATIONS 

A generalized conformal transformation is the transformation from the Lagrangian 
G(x,v) to the Lagrangian G(x,v) where 

G(x, v) = G(x, v) J(x, v)~ l . 

Here J(x, v) is a non-zero function, smooth and defined in a neighborhood of the 
cone G(x, v) = 0. We have the homogeneities, valid for any real t > 0: 

G(x,tv) = t k G(x,v) 

G(x,tv) = t p G(x,v) 

J(x,tv) = t q J(x,v) 

where 

p — q = k 7^ 1, p ^= 1. 

Denote by V the dynamical vector field with respect to the transformed La- 
grangian G. 
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Lemma 17. The dynamical vector field V transforms via V = V + bH + GT where 
b = —(p — l)~ 1 J~ 1 V(J) — (k — 1) J V( J) and T is some vertical vector field. 

Proof. Since both V and V satisfy Lemma[9j X(V — V) = 0, and therefore the two 
vector fields can only differ by a vertical vector field. But working modulo G = 0, 
since V are dynamical vector fields on Jif, V must be in the span of V and if, and 
so 

V = V + bH (mod G = 0) 

for some function b. 

It remains to show that b = -(p - 1) ^ J ^V (J) = -(k - l)~ 1 J~ 1 V(J). By 
definition of a, a — J a + G DJ, and so pulling back to Jff gives 

da = dJ A a + J da (mod G, dG). 

Contracting with V gives on the one hand V-ida = V(J)a since V-ia = and 
V Ada = — (k — l)dG = 0. On the other hand, using V = V — bH gives 

V-ida = —b(p — l)a = —b(p — l)Ja 

because V-ida = — (p — l)dG = and H_ida = (p — l)a. Combining these two 
calculations gives b — — (p — J~ 1 V( J), as claimed. By symmetry, —b = — (k — 
l^JViJ- 1 ) = (k - ly^J^ViJ), and so b = -(k - l)" 1 J" 1 ^ J" 1 ) as well. □ 

Lemma 18. Jz^A = £f v \ -bX- Db® H - a®T + G£? T X. 

Proof. 

^yX(X) = [V, XX] - X[V, XX] 

= [V + bH + GT, XX] - X[V + bH + GT, X] 

= ^fvX(X) + b£f H X(X) - (D x b)H + G^ T X{X) - (D X G)T + X(G)XT 

= (Jz^VA -bX- Db® H + GJ^tX -a®T) (X) 

where the last equality follows since S£rX = —A, DxG = a(X) by definition, and 
AT = since T is vertical. □ 

6.1. Weyl tensor. Recall that Ajg> is the subbundle of TJrf? consisting of vectors 
X that annihilate a: a(X) = 0. Then VJf C Ajp and also the dynamical vector 
field V is a section of Aj#>. 

Definition 5. The umbral bundle is the vector bundle E over defined as the 
quotient of Aj%> by the kernel of gh\h x xhj e ■ 

The umbral bundle is a rank n — 2 vector bundle over . It is so named 
because in 52 the pullback of E along sections of can be regarded as a space of 
infinitesimal screens onto which an object, placed into the null geodesic spray, will 
cast a shadow. This interpretation is due to Sachs |25j . 

The kernel of gh\K x xA^ is the subspace of spanned by and the dy- 

namical vector field V. Indeed, gh(V, X) = (k — l)a(X), which vanishes if X £ Aj#>. 
Thus ker((7ft|A^xAj^) contains V. It also contains VJff, since the kernel of the bi- 
linear form gh on the full tangent space TTM is VTM . From rank considerations, 
ghlAje-xAye has degree of degeneracy at most n, and so the kernel must in fact be 
equal to VJf © spanK Thus 
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If X G A.jg>, denote by [X] the equivalence class of X in E. The tensor gh 
descends to a non-degenerate metric on E, via 

g E ([X] > [Y])=g h (X,Y). 

When it is restricted to A je> x kj#> , S vanishes if either argument is in the kernel 
of gh (by Corollary [1}. Thus by restriction S defines a section of E* (g> E*. Define 
an endomorphism S| : E -> E by setting g E {S%X, Y) = S(X, Y) for all X,Y £ E. 

Definition 6. The Weyl tensor W G T ^ ® E*) is the trace-free part of the 
restriction of the tidal force tensor S(X, Y) to (X, Y) G E x E. That is, 

W(X, Y) = S(X, Y) - f-L- tr S E ) g E (X, Y) 



for (X,Y) GE x E. 

Notice that we are now working on exclusively, and so G = 0. 

Theorem 6. The Weyl tensor depends only on .ffl C TTM' , not on the choice of 
defining function G. 

In other words, the Weyl tensor is conformally invariant with respect to the class 
of conformal transformations described at the beginning of 



Proof. By Lemma EH S E [X] = \ A l P(JfyX)X 



for X G Ajff. So to prove that 



W is conformally invariant, it is sufficient to compute P(Jz?~ A) on E, and then to 
neglect terms that are proportional to A, since these will only modify the trace. We 
shall therefore compute A modulo terms involving G, dG, a, since these are zero 

on E, modulo H since PH = H which is in VTM' and so also zero in E in which 
the vertical space is quotiented, and modulo V since PV = PV (mod H) = 0. We 
treat each term of Lemma [TBI in turn: 

(JS^JBVA) X = [V, (JSJVA)Jf] - (^ V X)[V,X] 

= [V + bH + GT, (3f v \)X\ - {Sf v X)\V + bH + GT, X] 

= ^X(X) + b(^ H ^ v X)(X) + & GT (5? V X){X) (mod G, dG, a, H, V) 

= £>$X(X) + J? GT (J? V X)(X). 

by homogeneity of V and A. Now J£ G t(JZvX)(X) vanishes modulo G for X G 
kerdG. Hence 

(J^JSVA) X = & V X(X). 

The second term is 

(-%(6A)) X ee V(b)XX + b(£> 9 \)X 

ee V(b)XX + b(£> v X)X - b 2 XX. 

The third term is 

{^ v {Db ®H))X = -D x b ® V ee 0. 
The remaining terms are zero, since they involve a and G: 

J£ V {GJ£ T X) EE 
(j2^(a®T)) s 0. 
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Thus we have 

J2*A = JS^A - &(jSf v A) + (6 2 - V(6))A. 
We now compute PJz?~ A. By the transformation law for Jz? y A, 

P = P - ^(6A + Db® H + a®T - G^ T a). 

Note that since j£f„A = 2P— Id, PJz? v \ = P and A^fyA = —A. Among the remaining 
terms are those involving the a <g> T term of P contracted with a term of Jz?~ A. Of 
these, it follows from a o P = that .SfyA preserves the annihilator of a, and so 
the term (a ® T)(6JfyA) term vanishes when restricted to A^. Since a o A = 0, 
the (a (8) T)A terms vanish. Finally, the term involving (a <8> T)(Jz?yA) vanishes on 
Aj? — ker a since the kernel of a is Lie derived along V. 
So, applying P to Jz?~ A gives 

PJ^?A = Pif^A - bP + (b 2 - 9(b)) X - ^bX^X - ^b 2 X 

= P££%x + (\b 2 - a - hx^ 2 x - bP 

It is now sufficient to show that the last two terms cancel; that is: 

AJz^A = -2P. 

We have 

XJ? 2 X(X) = X([V, [V, XX}} - 2[V, X[V, X}} + X[V, [V, X}}) 

= X[V,[V,XX}]-2X[V,X[V,X}] 

= X[V, [V, XX}} + 2X[V,X]. 

If X is in the image of P, the first term vanishes because im P = ker A, leaving only 
the second term which is — 2X. If instead X G kcrP, then JzfyA = —A. So 

X[V, [V, XX]} + 2X[V, A] = X[V, -X + X[V, A]] - 2AA 

= -A[V, A] + X[V, -XX} - 2XX 
= 2XX - 2XX = 

as required. Thus, in summary 

PSfZ X = P££ 2 X + (\b 2 - V{bfj X. 

Since the term multiplying A only modifies the trace of S, this completes the proof. 

□ 

7. RAYCHAUDHURI-SACHS EQUATIONS 

We review the Raychaudhuri-Sachs equation of standard general relativity. Let 
M be a spacetime manifold of dimension n > 3, equipped with an indefinite metric, 
g of signature (p, q). Denote the Levi-Civita connection of g by V. 

Let k be a null vector field that is nowhere zero and satisfies the equation of an 
affincly parametrized geodesic Vfcfc = 0. The integral curves of k are null geodesies 
that foliate M: that is, they constitute a null geodesic congruence. Associated to 
the vector field k is a natural vector bundle K of dimension n — 2 with a metric 
of signature (p — l,q — 1) (so Euclidean in the case where g is Lorentzian). This 
bundle consists of the (n — 2)-plane elements (or "screens") onto which the infinitely 
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near curves of the congruence would cast the shadow of an object. This bundle was 
introduced in this way by Sachs |25j . A precise definition of this bundle is in section 
1721 

The Raychaudhuri-Sachs equation then governs the rate at which this shadow 
expands (or contracts) as the screens advance along a particular geodesic of the 
congruence. A principal ingredient in the derivation of the equation is the notion 
of the divergence of A:, of which there are potentially several candidates (that turn 
out to agree): 

• The Lie derivative of the volume element of M along k. 

• The Lie derivative of a natural volume element for the bundle K. 

• The trace of the endomorphism Vfc. 

7.1. Notation and conventions. The curvature tensor R e T M (/\ 2 T* M®End(TM)) 
is defined by the relation 

R(X, Y)Z = (V x V y - VyV X - V [X ,Y])Z. 

The Ricci tensor is given by 

Ric(X, Y) = tr(Z i-> R(X, Z)Y). 

The metric defines an isomorphism between the tangent and cotangent bundles of 
M: define g : TM -> T*M by 

g(X):Y^g(X,Y). 

This is a self-adjoint transformation (by the symmetry of g) that is invertible (by 
the non-degeneracy of g). The inverse g _1 : T*M — > TM defines a metric g~ x on 
T*M by 

g-\a,P)=p{g-\a)). 

The metric will be used to convert vectors into covectors systematically using the 
"musical isomorphism" : 

• If X is a vector, define X — g(X). 

• If a is a covector, define a" = g~ 1 {a). 

The volume element of M is a density on M that is defined on a collection of 
vectors V\,...,v n by 

\Cl(vi,...,v n )\ 2 = \det[g(vi,Vj)]ij = i t ... tn \. 

This is a section of the density bundle | A" T*M\. If an orientation is given on 
M, then it is possible to choose a representative volume form, denoted fi, for the 
density In a distinguished oriented local orthonormal basi^] of vector fields 
X\ , . . . , X n , 

n(x 1: ...,x n ) = \n{x 1 ,...,x n )\ = l, 

extended by multilinearity. 

Since a local one-parameter group of diffeomorphisms must preserve orientation, 
the Lie derivative of |f2| along any differentiable vector field is well-defined. The 
divergence of a differentiable vector field X is defined by 

(divJC)|n|=JSfjc|n|. 



For a metric of indefinite signature, an orthonormal basis is any basis such that g(Xi, Xj) = 
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Now, for the Raychaudhuri-Sachs equations, assume in addition that k is hy- 
persurface orthogonal. This is equivalent to the condition that the distribution 
k 1 - = (k b )° C TM of (n — l)-planes annihilating k b be integrable in the sense of 
Frobenius: k b A dk b = 0. 

7.2. The umbral bundle of the congruence. Let k 1 - denote the distribution of 
(n — l)-planes orthogonal to k. Thus, at a point x G M, 

k^ = {veT x M\g(k,v) = 0} = (k b x )°. 

Lemma 19. The distribution k 1 - is Lie-derived along k. That is, if v 6 Tm^^), 
then S£ k v G T M (k x ). 

Proof. If v is a section of k 1 - , then g(k, v) = 0. So 

= k(g(k, v)) = g{V k k, v) + g(k, V k v) = g{k, V k v) 
= g(k,^kv)+g(k,W v k) 

= g(k,& k v) + -v(g(k, k)) = g{k,^ k v) 

so J£ k v is also in k , as required. □ 

Note that k is a section of k 1 - , since it is null. Therefore the following definition 
makes sense: 

Definition 7. Let K be the quotient bundle K = fc^/spanfc. 

If a small object is placed in the path of the congruence k, then the bundle 
K naturally describes a family of screens onto which the shadow of an object is 
cast. Hence, this is the umbral bundle for the null geodesic congruence k. It is 
the pullback of the umbral bundle defined in §6.1l bv the section k of the null cone 
bundle Jf; see 33 for more details. 

Let [v] denote the equivalence class of v 6 k^ modulo k. Since k 1 - and spanfc are 
both Lie derived along k, the Lie derivative S£ k descends to a differential operator 
on the quotient K , by setting 

The Lie derivative extends to a unique derivation on the tensor algebra of K that 
commutes with tensor contraction. 

The metric g in TM induces a bilinear form g k ± on k , and the vector k is in 
the kernel of g k ± . Hence g k ± descends to a bilinear form on K via the rule 

9k{[X],[Y]) = g k x(X,Y). 

The bilinear form gx is a metric of signature (p — 1, q — 1) on K. 

The tidal force along k is the endomorphism S 1 " : TM — > TM given on vectors 
X by 

S^X = R{k,X)k. 

bmce S^k = and the image of is orthogonal to k, S 1 " induces an endomorphism 
of K via 

S* K [X] = [S*X]. 
The bilinear form S on TM and Sk on K given by 

S(X,Y) = g(S*X, Y), S K ([X], [Y]) = g K {& K [X\, [Y]) 
are both symmetric, by the symmetries of the Riemann tensor. 
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7.3. Divergence. 

Definition 8. Let X be a vector field. The divergence of X, denoted AwX, is 
defined by the equation 

{d.wX)\tt\=&xP\ 

For the vector field X, define the cndomorphism VX of TM by VX : Y i-> VyX. 
Lemma 20. The divergence of X is the trace ofVX G I\f(End(TM)) 

div X = tr VX 

Proof. Let v±, . . . ,v n be a local basis of smooth sections of TM, and let a 1 , . . . ,a n 
be the dual basis of T*M, defined by a l (vj) — Sj. Let ft be the local section of 
A n T*M representing obtained by declaring the basis Vi to be positively oriented. 
First note that if a is a one-form and Y a tangent vector, then 

= Yj (a A fi) = a{Y)n -aAYjJ] 

so 

(7) aAYjQ. = a(Y)Q. 

By Cartan's identities, 

Jz? x O = d(X j = a * A ( x J n ) 



V A (V„,I)jn + ^a'AljV„,l] 

■j i 

i 

= ^ ^(V„ i X)0 = tr(VX)0 



by ©. □ 

Fixing an orientation on TM equips the bundle K with an induced orientation, 
and the associated volume forms are related by 

fcj n = k b a n K 

where j is the interior product. The validity of this equation does not depend on 
the choice of coset representative of flu modulo the ideal generated by k b , and so 
defines £Ik uniquely as a section of A n ~ 2 K* . If no orientation on M is specified, 
then this only defines a density \VLk\ in the determinant bundle | A" -2 K*\. 

Definition 9. Define div^- k by 

{dxv K k)\Vt K \ =££ k \Sl K \. 

Lemma 21. div^- k = divfc 

Proof. Working locally with an orientation on M, we have Aij f2 = k A VLk- So 

ift(bfi) =J? k {kAVL K ) 
kj^f k Q = kA^ k VL K 
(div k)k_s Q — (divK k)k A Q,k- 

□ 
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The image of the endomorphism Vfc lies in fc- 1 -, since g(\7xk, fc) = hX{cj(k, fc)) = 
0. Furthermore, k lies in the kernel of Vfc, since k is an affinely parametrized 
geodesic vector field. Therefore, Vfc descends to an endomorphism Vfc|if : K — » K . 

Lemma 22. For any p — 1,2, ... , 

tr((Vfcf ) = tr((Vfck) p ) 

Proof. In general, if A is an endomorphism of a vector space V whose image lies 
in a subspace W, then tr A = tr (A\w)- Since Vfc is a linear operator whose image 
lies in k x , tr(Vfc) p = tr((Vfc)*V )■ But (Vfc) p |^ = (Vfc| fei ) p , and so tr(Vfc) p = 
tr(Vfc|fcj- ) p . Now, since spanfc lies in the kernel of Vfc \ k ± , tr(Wk\ k ±) p = tr(Vfc|if) p , 
as required. □ 

7.4. Rate of change of the divergence. The purpose of this section is to com- 
pute the rate of change of the divergence of fc. Let R(k, —)k denote the endomor- 
phism i?(fc, -)fc:Ih> R(k, X)k. Then: 

Lemma 23. V fe Vfc = -(Vfc) 2 + R(k, -)k 

Proof. For a vector field X, 

(V fc Vfc)pf) = V fc Vxfc - V Vfc xfc 

= R[k, X)k + V x V fe fc + V [KX ]k - V Vfc xfc 

= R(k,X)k- W Vxk k 

= [-(Vfc) 2 +R(k,-)k](X) 

□ 

Lemma 24. 

fc(div fc) = - tr[(Vfc) 2 ] + Ric(fc, fc) 

= -tr[(Vfc| K ) 2 ] +Ric(fc,fc) 

Proof. The first equation follows by taking a trace from the previous lemma. The 
second equation follows from tr(Vfc) 2 = tr(Vfc|i^) 2 . □ 

Lemma 25. tr5 tt = trS^ 

Proof. The image of 5 1 " lies in k 1 - and the kernel of S 1 " contains fc. Thus the lemma 
follows by the argument of Lemma [22j □ 

7.5. Invariant decomposition. Let 

Vfc|if = Alt(Vfclif) + Sym (Vfc|x) H l — (divfc)Idx 

n — 2 

be the decomposition of Vfc|if into its irreducible components for the action of 
0(p — l,q — 1): the alternating, symmetric trace-free, and trace parts. Here the 
metric qk is used to identify End(K) with K*®K* in order to define the symmetric 
and alternating parts. 

For the next theorem, introduce the following notation, standard in the relativity 
literature when n = 4: 

• 8 = div fc is called the expansion of the congruence fc in the relativity 
literature 

• a = Sym (Vfc|if) is the shear tensor 
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• p = Alt(Vfc|if) is the rotation tensor 
Theorem 7. 

k(0) = - trip 2 ) - tv(a 2 ) + Ricffc, k) 

n — 2 

= - tr(p 2 ) - tr(a 2 ) - J— + trS". 

n — 2 

Proof. This is a restatement of Lemma [24] under the decomposition 

Vk\ K = p + (T + ^—eid K . 

n — 2 

The absence of cross-terms owes to the orthogonality of the different irreducible 
representations of 0(p — 1 , q — 1). The second equality follows from the definition 
of Si ' ' □ 

7.6. Hypersurface orthogonality. If k is hypersurface orthogonal, then the dis- 
tribution k 1 - = (k b )° is integrable in the sense of Frobenius, and therefore dk b = 
(mod Ar). 

Lemma 26. If k is hypersurface orthogonal, then tr(p 2 ) = 0. 

Proof. If k is hypersurface orthogonal, then there exists locally a one-form \x such 
that dk b = p Ak b . Since A is a geodesic vector field, kjdk" — 0, and since k is also 
null k b (k) = 0, so p(k) — as well. Now 

tx{p 2 ) = rffc b (fc,/i») = 

as claimed. □ 

Theorem [7] becomes the Ray chaudhuri- Sachs equations: 

Corollary 2. If k is hypersurface orthogonal, then 

k(6) = -tr(o- 2 )--^+trS t . 

n — 2 

7.7. Raychaudhuri effect. 

Lemma 27. Suppose that J is a vector field that Lie commutes with k. Then J is 
a Jacobi field along any integral curve of k. 

Proof. Covariantly differentiating = [k, J] = \7 k J — Vjfc along k gives 

= V^J- V fc Vjfc 
= V 2 k J-R{k,J)k 

which is the Jacobi equation □ 

In particular, since k is hypersurface orthogonal, there are n — 2 (Jacobi) vector 
fields Ji, . . . , J n -2 that are orthgononal to k, Lie commute with k, and are linearly 
independent of k. On passing to the quotient, these Jacobi fields define a basis of 
K. Pick such a basis, and let \k = \&k(Ji, ■ ■ ■ , J n -2)\- 

In the Lorentzian case of a space-time of n-dimensions, the signature of the 
metric gE of the bundle E is either positive or negative definite, according as g has 
signature (n — 1, 1) or (l,n — 1). Thus in the Raychaudhuri-Sachs equations, the 



2(5 



JONATHAN HOLLAND AND GEORGE SPARLING 



trace tr(a 2 ) is non-negative, and it is zero if and only if a = 0. Thus Corollary [2] 
gives 

k(9) < tr5 tt , 

or equivalently, 

jSf fc 2 A K < tiSh K . 
The null positive energy condition is the condition 

Ric(n, n) < for all null vectors n. 
So when the null positive energy condition holds, 

^X K < 

Note that this equality only requires that Ric(fc, k) < be valid for the particular 
tangent vectors along the given null geodesic. 

Now suppose that 9 < at some point Xo of the congruence. By definition of 9, at 
that point Jzfk^K = OXk < 0. Then Xk will become zero along the geodesic tangent 
to k through xq at some time prior to the finite affine parameter t = —(ji — 2)/9(xq). 
Since k is hypersurface orthogonal, the vectors J\, . . . , J„_2 span the tangent space 
of this hypersurface up to the point where the volume Xk degenerates to zero. At 
or before that point, the geodesic in question must have a conjugate point. The 
existence of this conjugate point is the key to the proof by Sir Roger Penrose [21] 
of his singularity theorem. 

8. The geometric Raychaudhuri-Sachs theorem 

In this section, we lift the geometry underlying the Raychaudhuri-Sachs theorem 
to the bundle Jf? and at the same time generalize it to regular causal geometries. 
We first recall some basic sheaf theory. 

Let p : Y — > X be a (continuous) map of topological spaces. 

• A point y 6 Y is said to be a sheaf point if and only if there exists an open 
set U y in Y, such that y S U y and such that the restriction of p to U y is a 
homcomorphism onto its open image p(U y ) in X. 

• The sheaf space S p c Y of p is the collection of all its sheaf points, with 
the induced topology. Note that S p is an open subset of Y. 

• The triple (Y,X,p) is said to be a sheaf if and only if p is surjective and 
S P = Y. 

• The triple (Y, X,p) is said to be a stack if and only if p is surjective and S p 
is dense in Y, i.e. the closure S p — Y. 

• The triple (Y, X,p) is said to be a branched cover if and only if it is a stack 
and both Y and X are Hausdorff topological spaces. 

For example: 

• Put 8 1 = {(x, y) G R 2 : x 2 + y 2 = 1}, the unit circle in the plane. Then the 
map e : K — > 8 1 given by the formula e(t) = (cos(i), sin(t)), for any t e I 
makes (R, S^e) a sheaf. 

• Consider the complex parabola Y = {(x, y) G C 2 : y 2 — x} and let p(x, y) = 
x G C for any (x,y) G Y. Then the triple (Y, C,p) is a stack, with S p = 
Y — {(0, 0)} and is a branched cover. 
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If (Y, X,p) is a stack, and if U is an open subset of X, then a section s of the stack 
over U is a map s : U — > X, such that p o s = idu- 

The key concept we need is that of a Sachs manifold. Let M be a manifold of 
dimension n and let S? C SM be a regular causal geometry. Denote the natural 
surjection from Sf to M by p. Denote a representative one-form of the contact 
structure of Sf by a^. Put = and denote by a# a representative one- 

form of the contact structure of J$f . Also denote by q the natural surjection from 
Jf to M. 

• A Sachs manifold for the causal geometry is a smooth submanifold 5 of 
^ of dimension n, such that: 

— The triple (S, M,p\s) is a branched cover. 

— S is ruled by (unparametrized) null geodesies: i.e. the null geodesic 
spray V of Sf is everywhere tangent to 5. 

— S is hypersurface orthogonal: the restriction of the three-form u<$doi<$ 
to S vanishes identically. 

• An affine Sachs manifold for the causal geometry is a submanifold T of 
J^ 9 of dimension n, such that: 

— The triple (T, M,q\f) is a branched cover. 

— T is ruled by affinely parametrized null geodesies: i.e. the null geodesic 
spray V of 34? is everywhere tangent to T. 

— T is hypersurface orthogonal: the restriction of the three- form a,^pda 
to T vanishes identically. 

• A Sachs section for a causal geometry 34? over an open set U C M is a 
section of a given Sachs manifold, whose domain is U. 

• A Sachs congruence on an open subset U C M is the foliation of £/ by 
the null geodesies giving the foliation of a Sachs section. Note that the 
Sachs congruence is automatically hypersurface orthogonal, with normals 
the (null) tangent vectors to the foliation and the congruence and the sec- 
tion determine each other uniquely. 

In the special case of a standard space-time, the Sachs congruence exactly agrees 
with the congruence needed for the Raychaudhuri-Sachs equation and we see that 
in that case the affine Sachs manifold is simply the natural lift to the tangent bundle 
of the Sachs congruence, so we have a natural generalization. 

Now let a Sachs section s : M — > 3f be given. The connection P : T3V — > V34? 
defines an endomorphism P s 6 End(TM) given by 

Let k be the tangent vector field of the congruence, so s*fc = V. The tensor 
pulls back under s to a metric g s — s*gh on M. Moreover, k b — s*a. Since 
gh{V 1 V) — k(k — 1)G, it follows that g s (k, k) = since s is a section of 34? where 
G = 0. 

The bundle K is defined as before as fc^/spanfc, where k 1 - is the orthogonal 
complement of k with respect to the metric g s . This is naturally isomorphic to 
the pullback under s of the umbral bundle E defined in section 16.11 The metric 
g s induces a metric gx on K, which is of definite signature if g v has Lorentzian 
signature. Let V denote the Levi-Civita connection of g s - The Lie derivative J?k 
preserves ker k b , by Lemma IT9l Likewise the Lie derivative extends to all associated 
tensor bundles. 
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Lemma 28. P s = Vfc where V is the Levi-Civita connection associated with the 
metric g s . In particular k is an affinely parametrized geodesic with respect to the 
connection V. Moreover, the pullback of the tidal force tensor along s is the sectional 
curvature o/V in the direction of k: 

S(s*X,s*Y) = g s (R(k,X)k,Y) 

where R is the Riemann tensor associated to V. 

Proof. The proof of the first claim proceeds by verifying that the two tensors have 
the same skew and symmetric parts. On the one hand, 

(JSfkSaX*, Y) = k(g s (X, Y)) - g s {V k X - V x k, Y) - g s (X, V k Y - Vyfc) 

= g s {W k X, Y) + g s {X, V k Y) - g s {V k X - V x k, Y) - g s (X, V k Y - Vy 
= g s (V x k,Y)+g s (X, Vyfc). 

On the other hand, 

(& k g.)(X,Y) = (^ v g h )(s^X,s t Y) = g v (Ps*X, Xs,Y) + g v {\s*X, Ps*Y) 

= ghiX^Ps.X^.Y) + g v (s*X,\~ 1 Ps*Y) 
= g s (P s X, Y)+g s (X,P s Y). 

This shows that Vfc and P s have the same symmetric part. 
For the skew part, on the one hand 

2(s*da)(X,Y) = 2dk\x,Y) 

= g s (V x k,Y)-g s (X, Vyfc) 

and on the other hand 

2(s*da){X,Y) = 2da(s*X,s*Y) = g v (Ps*X, Xs,Y) - g v (\s,X, Ps*Y) 

= ghiT 1 Ps*X, s*Y) - ghis.XX 1 Ps*Y) 
= g s (P s X,Y)-g s (X,P s Y). 

Since P s = Vfc, V^fc = P s k = A 1 PV = since V is horizontal for the Ehres- 
mann connection P. Hence fc is an affinely parametrized geodesic. 
For the final claim, Theorem |4] implies that it is sufficient to prove 

i(j^ 2 5s )(X, Y) = g s (V x k, Vyfc) - g s (R(k, X)k, Y) 
since P s = Vfc by the first part of the lemma. The identity 

Jzf fe Vfc = V fc Vfc 

holds, so 

l -{^g s ){X, Y) = g s (V x k, Vyfc) + ^g s ((V k Vk)(X), Y) + ^g s (X, (V fc Vfc)(F))+ 

+ ^ s (v Vx ^,y) + -g s (X, V Vrfc fc) 

- g s (V x k, Vyfc) + i (g s (R(k, X)k, Y) + g s (X, R(k, Y)k)) 
= g s {V x k, Vyfc) + g s (R(k, X)k, Y) 
by the symmetries of the Riemann tensor. □ 
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The operator S\, : ~k^} m ,TM — > ■k' t 1 mi TM defined in £ 16.11 when restricted to the 
section s defines an operator S* : TM — > TM. By the previous lemma, S\{X) = 
R s (k, X)k. Moreover, as in fJ7J the image of S| lies in k 1 - and its kernel contains fc, 
so 5| descends to and operator : K — > K. Moreover, tr S\ = tr S\ { = Ric(fc, k) 

As in £ 17.31 the divergence of k can be defined in several equivalent ways. If |0| 
is the canonical density associated to the metric g a , then 

jS? fe |fi| = (divfc)|0|. 

If \£Ie\ is the canonical section of the determinant line bundle | A" -2 E\, then 

Sf k \Q E \ = (div B fc)|fi B |. 

Alternatively, the divergence can be defined as the trace of Vfc = P s , or the trace 
of Vk\ E = P s \e- The results of §7.31 imply that these are equal: 

Lemma 29. 9 = divA; = div B A; = tr(Vfc) = tr(P s ) = tr(V fe | B ) = tr(P s \ E ) 

The proof of Theorem [7] goes through as in fJT] 

Theorem 8. Let 

P s \ E = Alt P s | B + Sym P s \ E + Id ^ 
9 

= p + a H Mb 

n — 2 

6e t/ie decomposition of P s into its irreducible 0(p — l,q — 1) components. Then 

k{9) = - tr(p 2 ) - tr(a 2 ) - + tr S*. 

n — 2 

8.1. The Lorentzian case: the geometric Raychaudhuri Sachs effect. Now 

consider the case that the fibre metric g v in VTM' is Lorentzian, which implies in 
turn that g s is also Lorentzian, and so the metric of the bundle K has positive 
or negative definite signature. Then the quantity tr(<r 2 ) of the Raychaudhuri-Sachs 
equation is non-negative. Also impose the positive energy condition: tr 5" < 0. As 
in $71 let J%, ... , J n -2 be a collection of vector fields orthgonal to k that commute 
with fc, and set \k — |Oif(Ji, . . . , J n -2)\- Then 

^k^K — GXki ^k^K < 0. 

Now if at a point of the congruence we have 9 < 0, then it follows that the graph of 
\k along the (affinely parametrized) null geodesic through the point is decreasing 
and concave down, so Xk reaches zero in finite affine parameter time in the future. 
So we have the theorem: 

Theorem 9. Let X be a given null geodesic in M that is future complete, so its 
affine parameter ranges to positive infinity. Suppose that everywhere along X the 
positive energy condition tr Sfl < holds. Suppose there is a section of a Sachs man- 
ifold, defined in a neighborhood of X , such that X is a member of the congruence 
foliating the Sachs manifold. Then the divergence of the congruence is everywhere 
non-negative along X . 
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Appendix A. Notational conventions 

• If M is a manifold and E — )• M is a bundle, then the projection is denoted 
by tte ■ The space of smooth sections is denoted by Tm{E). 

• If A is a vector space and S C A* is a subset, then the annihilator of S, 
denoted by S° , is defined by 

S* = {x e A\a(x) = for all a e 5}. 

• Arbitrary vector fields are denoted by uppercase latin letters at the end of 
the alphabet: W, X, Y, Z. In Section [71 k is used to denote a null geodesic 
vector field. Throughout the paper, V denotes the null geodesic spray and 
H the homogeneity operator. 

Appendix B. Coordinate calculations 

B.l. The Hamiltonian-Lagrangian approach to the dynamics. Let Jt? have 
the local defining equation G(x, v) — 0. Here the smooth function G(x, v) is denned 
over an open set U of TM , which is invariant under scaling: i.e. (x, v) G U implies 
that (x, tv) € U, for any positive real t. Further we may take the function G{x, v) 
to be homogeneous: G(x,tv) — t k G(x,v), for some real k and any t > and any 
(x, v) in the domain of G(x, v). For convenience, we henceforth assume that k =/= I. 

In local co-ordinates (x a ,v a ), where a = 1,2, ... ,n, denote by d a and D a , the 
derivative operators: 

d a = — D a = — 

a dx a ' dv a 

These operators mutually commute. Define the following quantities: 

G a = d a G, g a = D a G, g a b — D a D G — gh a , g a bc = D a DbD c G = g( a bc)- 

Note that we have: 

v a g a = v a D a G = kG, v a g ab = (k - l)g b , v a g ab c = {k - 2)g bc . 

Then, by the last section, a representative contact one-form for the dynamics is the 
following one-form a, with exterior derivative (3, considered on the space G = 0: 

a = g a dx a , (3 = da = -(D{ a G b ])dx a dx b + g ab dv a dx b . 

For V = V a d a + U a D a to be a dynamical vector field, we need the conditions: 

= G, = dG = G a dx a + g a dv a , = V(G) = V a G a + U a g a , 

Q = V.a = V a g a , 

= V.P- t(k - l)a = -V a g ab dv b + (U a g ab - 2V a D [a G b] - t(k - l)g b )dx h . 
Here t is a scalar function. We infer the relations, for some scalar function s: 

V a g ab - s(k - l)g b , 

U a g a b - 2V a D [a G b] - t(k - 1)56 = s(k - l)G b . 
Henceforth we assume that g a b is invertible (the regular case). Then we have the 
general solutions: 

(V a ,U a ) =t{0,v a ) + s(v a ,u a ), 
u a g ab = 2v a D [a G b] - (k - l)G b = G b - v a d a9b . 
So a basis for the collection of dynamical vector fields is the pair {H, V}: 

H = v a D a , V = v a d a + u a D a , u a g ab = G b - v a d a9b . 
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Note that we have the required relation: 

V(G) = v a G a + u a g a = v a G a + (k - l)- 1 u a g aA v b 

= (k- l)-\(k - l)v a G a + (G b - v a d a9b )v b ) 

= (k-l)- 1 (kv a G a -v a d a (v b g b ))=0. 

Note that H is the homogeneity operator in v a . Also u a is homogeneous of degree 
2 in v a . Passing down to the sphere bundle §M, only the direction field V survives, 
giving a foliation of 'S by the dynamical curves, the (maximally extended) trajec- 
tories of V. We call these dynamical curves null geodesies. Note that the space 
of null geodesies, denoted jY , which has dimension 2n — 3, itself carries a natural 
contact structure induced by that of . 

Summarizing, the hypersurface uniquely determines its dynamics, which gives 
a foliation of Sf by a 2n — 3 parameter set of null geodesic curves. These curves are 
precisely those that annihilate the contact structure ay and preserve the contact 
structure (up to scale). 

B.2. The Ehresmann connection and its characterization. Introduce the 
quantity U b a and the horizontal vector fields H a : 

U b a = -^D b u a , H a = d a -U a b D b . 
Note that U b a is homogeneous of degree one in v a . Also we have the relations: 

v b U b a = -u a , 
V(g a ) = v b d b g a + u b g ab = G a , 
U b a g a = l -{-D b {u a g a ) + G b - v a d a g b ) = l -{D b {v a G a ) + G b - v a d a g b ) 

= \{v a (D b G a -d a9b ) + 2G b ) = G b . 

This last relation shows that the vector fields H a are intrinsic to the dynamical 
surface G = 0: 

H a (G) = (d a - U b D b )G = G a - U b g b = 0. 
Note that we have the simple formula determining the null geodesic spray V in 
terms of the horizontal vector fields H a : 

V = v a H a . 

The intrinsic tangent vector fields of G — are spanned by {H a , L bc }, where L bc = 
g[ b D c ^ . Dually, wc introduce the basis of one-forms: 

9 a = dx a , 4> a = dv a + U b a dx b . 
These are subject to the single linear relation, valid in the space G = 0: 
g a r = g a dv a + g a U b a dx b - g a dv a + G a dx a = dG = 0. 

For Sf we may use the forms 9 a and 4> ab — v^ a (j) b ^ as a spanning set. 

The horizontal vector fields H a form an Ehresmann connection over M, for the 
space Sf . We now show how to characterize this connection and the corresponding 
curvature, uniquely. 

The curvature generator is (symmetric tensor product here): 

Go = \ 2 g ab dx a dx b . 
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We take the Lie derivative of Go along V, giving the symmetric tensor Q\\ 

Qi = J^vQa = g a bdv a dx b + ^V(g ab )dx a dx b . 
Next we Lie derive Q\ along V, giving the symmetric tensor Q 2 : 

Qi = &vGi = 2V{g ab )dv a dx b + ^V 2 (g ab )dx a dx b + g ab dv a dv b + g ab dx b du a 

= 2V(g ab )dv a dx b + ^V 2 (g ab )dx a dx b + g ab dv a dv b + g cb dx b dv a D a u c + g bc dx a dx b d a u c 
^ g ab dv a dv b + dv a dx b (2V(g ab ) + D a (u c g bc )-u c g abc ) + ^dx a dx b (V 2 (g ab ) + 2g bc d a u c ) 
= g ab dv a dv b +dv a dx b {2v c d c g ab +u c g abc +D a (d b G-v e d e g b )) + ^dx a d 

= g ab dv a dv b +dv a dx b (2v c d c g ab +u c g abc +d b g a -d a g b -v e d e g ab )+^dx a dx b (V 2 (g ab )+2g bc d a u c ) 

= g ab dvUv b +dv a dx b (V(g ab )+d b9a -d a g b )+^dxUx b (V 2 (g ab )+2d a (u c g bc )-2u c d a g b ^ 

= g ab dv a dv b +dv a dx b {V{g ab )+d b g a -d a g b )+^dx a dx b {V 2 {g ab )+2d a d b G^ 

^g ab dv a dv b +dv a dx b (V(g ab )+d b g a -d a g b ) + ^dx a dx b (V 2 (g ab )-2V(d a g b )+2d a d b G). 
Now recall that: 

u a g ab = d b G - v e d e9b , U b = -2- 1 D a u b , 
v a U b = -u b , U a b g b = G a , 
4> a = dv a + U b a dx b , 9 a = dx a . 

Then we have: 

Uab = U a c g bc = -2- 1 D a (u b ff 6c ) + 2- 1 u b g abc 

= -2- l D a (d b G - v e d e g b ) + 2- 1 u b g abc 

= -2- x d b g a + 2~ 1 d a g b + 2- 1 v e d e g ab + 2" V<7 afec 

= 2- 1 (V(g ab ) + d a g b - d b g a ), 

v a U ab = -u b , U ab v b = G a . 

Substituting into our expression for Q2 and completing the square on the terms 
involving dv a , we have: 

G2 = gabdv a dv b +dv a dx b (V(g ab )+d b g a -d a g b )+^dx a dx b (V 2 (g ab )-2V(d a g b )+2d a d b G) 

= g ab dv a dv b + 2dv a dx b U ba + ^dx a dx b (V 2 {g ab ) - 2V(d a g b ) + 2d a d b G) 

= g ab (dv a +U c a dx c ){dv b +dx d Uj>)+^dx a dx b {V 2 (g ab )-2V(d a g^ 

= g ab r^ b + \dx a dx b {V 2 {g ab ) - 2V(d a g b ) + 2d a 8 b G - 2g cd U a c U b d ). 
So we have now: 

Q2 = g ab r^ b + o a e b s abl 

S ab = 2- 1 (V 2 (g ab ) - 2V(d a g b ) + 2d a d b G - 2g cd U a c U b d ) = S ba . 
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The tensor S a b generalizes the sectional curvature tensor of Lorentzian geometry. 
Back substituting into Q\, we have, since the skew terms in U ba cancel with the 
symmetric dx a dx b : 

Qi = g ab o a <p b . 

Now suppose we modify the Ehresmann connection <f) a to ip a — cf) a + 9 b u bc g ac , for 
some u ba - Then we have Q\ — g a bi ,a Q h provided that u ab = u ba . Also expressed in 
terms of ifi a and 9 a , we have: 

G 2 = 9abV^ b - ^ a e b u ba + e a e b u bc u bd9 cd + e a e b s ab . 

We see that the cross term in tp a 9 b vanishes if and only u ab — 0. We have proved 
that the Ehresmann connection is uniquely characterized by the absence of 6 a (j> b 
terms in Q2, so by the decomposition: 

G 2 = + o a e b s ab . 

Then, as a bonus, we have the formula: 

Gi = gab8 a 4> h - 

B.3. The curvature and its relation to the sectional curvature. For the 

curvature of the Ehresmann connection, we have: 

[H a ,H b ] - [d a - U a c D c , 8 b - U b d D d ] 
= -2(d [a U b f - U { d D w U b ^)D c = -2R ab c D c , 
R ab c = d [a U b] c U [a d D w U b f = H [a U b f. 

Note the relation: 

R ab c g c = d [a (U b] c g c ) + U [a c d b]9c ~ U [a d D w (U b] c g c ) 

= d [a G b] + U [a c (d b] g c - D\ c \G b] ) = 0. 

We extend the vector fields H a and D a to act on forms, as derivations of degree zero, 
by requiring that they annihilate the forms 9 a and 4> a . Introduce the derivations of 
degree minus one, denoted S a and e a , dual to 9 a and cf) a , which obey, in particular, 
the relations: 

s a e b - e a 4> b - s b , s a <p b = e a e b = o. 

The operators S a and e a mutually anti-commute and commute with H a and D a . 
Since d9 a = 0, we have the expression for the exterior derivative: 

d = 9 a H a + 4> a D a + {d^ a )e a . 

Define the covariant exterior derivative: 

d = dx a H a . 

Note the relation: 

dv a = dx b H b v a = -U a , U a = dx b U b a . 

Then, using forms, we have the curvature two-form operator R and the curvature 
two- form R c : 

R=-d 2 = -{dx a H a f = (dU c )D c = R C D C , 
R c = dU c = dx a dx b R ab c = dx a dx b (d a U b c - U d D d U b c ) = dx a dx b H [a U b] c . 
The Bianchi identity for the covariant derivative, d, is: 

dR b = d 2 U b = -R a D a U b . 
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Written out this is: 

dx a dx b dx c {H a R bc d ) = -dx a dx b dx c R ab e D e lI c d , 

H [aRbc] d + R [ab D c}U e d = 0. 

Introduce also the Lorentz generators, which are tangent to 'S and which, to- 
gether with H a span the tangent space to ( S: 

L ab = ^(9aD b - gbDa) = 9[ a D b] . 

We have the commutators: 

[Lab, L cd ] = [g[aL>b],9[cL>d]} = 9[agb\[cL>d] ~ 9[c9d][aL>b] — 9[c[bga]D d ] ~ CJ{a[d9c}L>b} 

= -^9[c[aLb\d\ = ^9{a[ C Ld]b}- 

[H a ,L bc ] = ^[H a ,g b D c - g c D b ] = (H a g [b )D c] + g [b (D c] U a e )D e 

= (H a g [b )D c] - (D [b U^)g c] D e 
= (H a9[b )D c] - 2(D [b U^)L c]e - g e (D [b U^)D c] 
= (H a g [b )D c] - 2(D [b U^)L c]e - (D [b (g e U^))D c] + U a[b D c] 
= (H a g [b )D c] - 2(D [b U ]a f)L c]e - (D [b G H )D c] + U a[b D c] 
= (H a9[b )D c] - 2(D [b U^)L c]e - (d a g [b )D c] + U a[b D c] 
= -U e a g e[b D c] - 2(£> [b [/ |a f)L c]e + U a[b D c] = -(D [b U^)L c]e . 

So we have: 

[H a ,L bc ] = -2(D [b U\ a \)L c]e = -2(D a U [b e )L c]e = U a[b L c]e , 

U ab c = -2D a U b c = D a D b u c = U ba c . 
Our second version of the sectional curvature is defined as follows: 
T b = v a 6 a R b = 2v c dx a R ca b = v a 6 a dU b 
= v a [8 a ,d]U b - v a 8U b = v a H a U b - v a dU b 
= v a H a U b - d(v a U a b ) + U b dv a 
= VU b + du b - U a U b = dx a T b , 
T b = 2v c R ca b = V{U b ) + H a u b - U a c U b . 
Lowering the upper index of T a b , we have: 

Tab = T a c g bc = V(U a c )g bc + g b cH a u c - U a c U cb 
= V(U ab ) - U a c V(g bc ) + H a {u c g ch ) - u c H a g bc - U a c U cb . 
So now we compute: 

2T ah - 2S ab = 2T ab - V 2 (g ab ) + 2V(d {a g b )) - 2d a dbG + 2g cd U a c U b d 
= 2V(U ab ) - 2U a c V(g cb ) + 2H a {dbG - v c d c g b ) - 2u c H a g cb 
-2U a c (U cb - Ubc) - 2d a d b G - V 2 (g a b) + 2Vd (a g b ) 
= V 2 (g ab ) + V(d a g b - d b g a ) - 2U a c V(g cb ) + 2H a (d b G - v c d c g b ) - 2u c H a g cb 
-2U a c (U cb - U bc ) - 2d a d b G - V 2 {g ab ) + V(d a g b + d b g a ) = 2Y ab , 
Y ab ^V{dagb)-U a c V(g cb )Hda-U a e D e )(d b G-v c d c gb)-u%d a -U a e D e )g cb -U a c (d c g^ 
= V{d a gb)-U a c V{g c b)-v c d c d a gb-U a e D e {dbG-v c d c gb)-u^da-U a e D e )g cb -U:{d c gb-dbg 
^u & D e dagb~U a c V{g cb )-U a e d b g e +U a e D e {v c d c gb)-u c {d a -U a e D e )g cb -U a c {d c g b -d b g c ) 



CAUSAL GEOMETRIES, NULL GEODESICS, AND GRAVITY 



35 



= u e dageb-U a c V(g cb )~U a e d b g e +U a e degb+U a e v c d c g eb ~u c (da-U a e D e )g cb -U a c (d c g b -d b g c ) 
= u e d a9eb - U a c V(g cb ) + U a e v c d c9eb - u c (d a - U a e D e )g cb 
= u e d a g eb - U a c u e D e9cb - u c (d a - U a e D e )g cb 
= u e d a g eb - u c d a g cb = 0. 
In going from the penultimate to the last line we used that the tensor D e g cb — g ecb — 
D e D c D b G is totally symmetric. So we have proved that the sectional curvature is 
simply the contraction with v a with the curvature, up to a constant factor: 

v a S a R» = dx a S b , S b = S ac g hc , 

^ va R ab = S b = S ab g ac . 
Note the relation implied by this, using the fact that R ab c is skew in a and b: 

v a S ab = 0. 

We check this relation directly: 

v a S ab = v a (V 2 (g ab ) - 2V{d {a9b) ) + 2d a d b G - 2g cd U a c U b d ) 
= v a V 2 (g ab ) - 2v a V(d (a9b) ) + 2v a d a d b G + 2g cd u c U b d 
^[v a ,u e D e ]V{g ab )+V{v a Vg ab )-2[v\u e D e ]{d (a g b) )-2V{v a d^ 
= -u a V(g ab )+V([v a ,V}g ab )+(k-l)V 2 g b +2u a d {a g b) -V(v a d a g b )-kVd b G+2v^ 

= -2 U a F( 5Q6 )- 5a ^( u a )-^G+2M a a (Q5b) - w Q a a ^-« a F(a a5t )+2« a a a a 6 G+2 5cdU c ^ d 

= -2u a V(g ab )-g ab V(u a )-Vd b G+u a d b g a ~v a V(d a g b )+2v a d a d b G+u a (V(g ab )+d^^ 
= -V(u a g ab ) - v a V(d a g b ) + v a d a d b G + u a (d b g a - d a9b ) 
= -V(8 b G - v c d c g b ) - v a V(d a9b ) + v a d a d b G + u a (d b g a - d a g b ) 
= -u a d b9a + (V(v a ))d a g b + u a {d b9a - d a9b ) = 0. 
Finally we wish to show that the sectional curvature S b c — g ac S ab determines the 
full curvature tensor. We start with the formula, just proved above: 

S b c - 2v a R ab c , 

Take the curl of both sides with D a giving: 

D [a S b] c = 2R ab c - 2v e D [a R b]e c . 

We need to analyze the last term, so we first recall the formula for 2R be c : 

2R be c = 8 b U e c - d e U b c - U b d D d U e c + U e d D d U b c . 

Then we have, using repeatedly the relation D a U b c = D b U a c and the fact that U a b 
is homogeneous of degree one in v a , so v c D c U a b = U a b : 

2v e D [a R b]e c = v e D [a d b] U e c + v e U [a 9 D b] D g U e c + v e {D [a U^)D\ g \U b] c 
= -v e D e d [a U b] c + U [a 9 v e D\ el D b] U g c + (v e D e U [a 9 )D\ g \UJ 
= -d [a U b f + U [a 9 D lgl U b f = -R ab c . 
So we have the required relation: 

D[ a S b f = 2R ab - 2v e D [a R b]e c = 2R ab - {-R ab ) = 3i? ab c , D [a R bc] d = 0. 

So we have proved that the sectional curvature determines the full curvature and 
conversely! Consequently we have the relations, using the Lorentz generators, in- 
stead of D a : 

L [ ab S c f = g[ a D b S c d = 6g [a R bc] d , L [ab R cd ] = 0. 
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B.4. Generalized conformal transformations. A generalized conformal trans- 
formation is the transformation from the Lagrangian G(x, v) to the Lagrangian 
H(x,v), where we have: 

G(x,v) — H(x,v)J(x,v)~ 1 . 

Here J(x, v) is a non-zero function, smooth and defined in a neighbourhood of the 
cone G(x, v) = 0. We have the homogeneities, valid for any real t > 0: 

G{x,tv) = t k G{x,v), H(x,tv) = t p H(x,v), J(x,tv) =t q J(x,v), 

p-q = k^ 1, p^ 1. 

Write J = e s , where v a D a s = q. Taking derivatives, we have: 

H = e s G, 

hi = e s (g l + SiG), 

ji = c Sj, jij = e (sij -\- Sj5j), 

% = e s (( s ij + SiSj)G + Sigj + s J g l + gij). 

Here we have: 

[dii Hi hi, Si] — Di[g, h,j, s], [yijiJiji hij, Sij] — DiDj[g, h,j, s] — [gji,jji, hj%, Sji\. 
We also put: 

[Gi, Ji, Hi, Si] — di[G, H, J, S], Ji = e s Si. 

Note that we have: 

//, = JjG + JGj = e s (Gj + SjG). 
We assume henceforth that the matrix gij is invertible. Using the relations 
v% 9ii — — l)ffj> v ^9j — kG and v l v 3 — (k — l)tPgj — k(k — 1)G, we have: 

e~ s hij = (s^ + SiSj)G + s t gj + s J g l + g tj 

= g Pg (Sf + {k- 1)- V«i)(<J? + {k- ly^Sj) + G{ Sij - (k - l^SiSj). 
On the cone G = 0, we deduce the relations: 

e-'ha = g Pq MfM], 

M? = 5j + (k- 1)- Vaj, (Af- 1 )*. = 5) - (p - lJ-Va,-, 

On the cone G = 0, the matrix e~ s hij is invertible, with inverse e 8 h li given by: 

e-fctf = ff w(M- 1 )*(M- 1 )j. 

Then, on G = 0, the signatures of <7jj and hij are equal if J > and opposite, if 
J < 0. Note that it follows that hij is invertible in a neighbourhood of G = 0, with 
inverse still denoted h % K Contracting the equation for e~ s hij with u l , using the 
fact that u a g a + v a G a = V{G) = 0, we get: 

hijU 1 = e s u l {s lJ +s l Sj)G + u\sie s gj + SjU l e s g l + e s dj(e~ s H)-e s v m d m {e~ s hj-s :j G) 
v m d m h J +e s u\sij+s l Sj)GWs i e s g^+s J u l e s gi-e s GS J +e% m d m (s J G)+v m S m e s {g J +SjG) 
= Hj - v m d m hj + V(s)e s gj + e s G(u l Sij + ijsiSj - Sj + v m d m Sj + v m S m Sj) 
= Hj - v m d m hj + V{s)e s gj + e s G(-Sj + V( Sj ) + V(s)sj) 
= Hj - v m d m hj + V{s)hj + e s G{-S 3 + V(sjj). 

So we have: 

u l = h ij (Hj - v k d k hj) + V(s)( P - 1)~V + e s Gh ij (V{ Sj ) - Sj). 
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We solve this equation for u a , order by order in G. To order G we see that: 

u a = 1 a + 5v a + Ge a , 
7° = h ab (H b - v k d k h b ), 5 = W(s)(k - 
V = W + 5v a D a + Ge a D a , 
W = v a 8 a + 1 a D a . 
Substituting, we need the relation: 
= h ij (H j -v k d k h 3 )+V(s)(p-l)- 1 v i +e s Gh ij (V(s j )-S j )-h ij (H j -v k d k h b j)-v i ^ 

= (W(s) + 5q + Ge a s a )(p - l)~ 1 v i + G(-e l + e s h ij (V( Sj ) - S 3 ) - v l W(s){k - l)" 1 . 
The terms not proportional to G cancel, so factoring out G, we get the relation: 
e l = e a s a (p - 1)-V + e s h ij (V(s 3 ) - S 3 ), 

e l = e a s a (p - 1)~ V + e s h ij (W( Sj ) - 5s 3 + Ge a s a3 - S 3 ) 
Next we have: 

e a = ( a + 7]v a + G9 a , 
^ = e sf l i i(W(s j )-6s j -S j ), 
V =(k-l)- 1 Cs a . 

= (C a +vv a + G8 a )s a (p-l)- 1 v l + e s h l i {W(s 3 )- 5s 3 + Ge a s aj -S 3 )~C- f]v l - G6 l 
= ((p - 1)77 + G6 a s a )(p - l)~ 1 v i + Ge s We a s aj - rtf - G9\ 
e i = e a s a (p-l)- 1 v i + e s h ij e a s aj , 
6 i = 6 a s a (p-l)- 1 v l + Ge s h ij 9 a s aj + e s h l3 {C + riv a )s a3 . 
Note that the residual term 9 % obeys the linear matrix equation: 

9 i - 9 a s a {p - 1)"V - Ge s h^9 a s aj = e s h^(( a + r]v a )s aj , 

(5 b a - s a (p - 1)- V + Ge 8 h bc s ac )9 b = e s h ij (C + r)v a )s aj , 
P{9) = Q, 
Pa =5 b a - s a {p - 1)- V + Ge s h bc s ac , 
Q a = e s h ij {C +rjv a )s a3 . 

When G = 0, P b has inverse the matrix M b given above, so the matrix P is always 
invertible in a neighbourhood of the surface G = and then 9 is given by 6 = P~ 1 Q. 
In particular the required solution exists and is unique in a neighbourhood of G = 0, 
as required. Also on the surface G = 0, 9 is given explicitly by the equation 
9 = M(Q) (i.e. 9 b = M a b Q a ). 

Summarizing, we have now the expansion: 

u a = 7 a + 5v a + G( a + Gr]v a + G 2 9 a 

Then on G = 0, we get: 

u a = 1 a + 5v a , 
U b a = -2~\D bl a + v a D b S + 55 a b + 9bC + v a 9 b v)- 

Here we have: 

7 a = h ab {H b -v k d k h b ), 
5 = W(s)(k-l)-\ 
C = e s hV(W(s 3 )-S S3 -S 3 ), 
V =(k-l)- 1 Cs a . 
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We apply these formulas to determine the conformal transformation of the sectional 
curvature. We have, evaluating on G = 0: 

T b = V{U b ) + H a {u b )-U a c U b . 

Using primes to denote the corresponding quantities computed for the conformally 
rescaled function H, we have: 

(T') a b = V'((U') b )+H' a {{u') b ) (U') a c (U') b . 
Now, on G = 0, we have: 

V = v a d a + u a D a = W + Sv a D a , 
V = W = v a 8 a + 1 a D a , 
(u' -u) a = -Sv a , 
H' a -H a = -(U'-U) b D b , 

(V - U) b = \{v b D a {5) + g a (( b + r,v b ) + SS b ). 
Using these relations, we have: 
(T') a b T b = W((U') b ) -(W + Sv c D c )U a b + H' a ((u') b ) H a u b (U') a c (U') b + U a c U c » 
= W((U' U) b ) 5U b + H' a ((u' u) b ) + (H' a H a )u b (U') a c (U') b + U a c U b 
= W((U' U) b ) 5U b + H' a ((u' u) b ) (V U) a c D c u b (U') a c (U') b + U a c U b 
= W((U' U) b ) 5U b H' a (Sv b ) + 2(U' U) a c U c b (U') a c (U') b + U a c U c b 
= W((U' U) b ) 5U b H' a (6v b ) (U' U) a c (U' U) b [V U, U] b a 
= W((U' - U) b ) + 5(U' U) b v b H' a (S) (V - U):(U> U) b [V U, U] b a . 
We write this expression out, working modulo multiples of v b and g a and working 
on G = 0: 

4(1* - T) b = 4W((U' U) b ) + AS(U' U) b A(U' U) a c (U> U) b A[U' -U,U] b a 
, 6 i) a («)+ 5a (C^ J ,^)+Mj)+# 2 ^(, c i}„(i) + ^)( Jc (^^))-2 9c C l [/ a c +2[/ c t , c ]} a (i) 
= 2{{V-5v c D c ){v b )-u b )D a {5)+2Q b {{V-5v c D c )(g a )-U a c g c )+{2{^ 

= (2V(S) + S 2 )S b . 
So, passing to the shadow space, we get the formula: 

(T')i = Ti + \(2V(6) + 6 2 )5i. 

We decompose and (T')£ as: 

{T't = {W't + \x'5l {WX=0, X ' = ^i T X, 

tp = wH + \x&Z, W% = 0, x = -±_t«. 

Then we have the key results: 

{w't = w£, 

X' = X + 2V(5) + <5 2 , 
(k - 1)6 = W{s) = V{s) -q8, 5 = (p- l)"V(s), 

x' = x + 2( P - iy V 2 ( S ) + (p - ir 2 (v») 2 

= X + 2(p - i)- V( J-V(J)) + (p - i)- 2 J- 2 (V( J)) 2 , 
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X' = X + 2(p - I)" 1 J-V 2 (J) - (p - ir 2 (2p - 3) J- 2 (F(J)) 2 . 

We have proved, in particular, that the trace-free part of the sectional curvature 
on the shadow space is conformally invariant. By definition, this conformally 
invariant part is called the Weyl sectional curvature and is a straight-forward gen- 
eralization of the sectional Weyl curvature of general relativity. Note that by judi- 
ciously choosing the conformal factor J, we may always assume (locally at least) 
that the trace of the sectional curvature is zero, so that T@ is trace-free. The 
residual freedom in conformal rescaling then obeys the relation 2V(S) + S 2 = 0, 
a first-order equation in <5 but second-order in the conformal factor. Finally note 
that in three dimensions, the sectional Weyl curvature vanishes identically, just as 
it does in three-dimensional metric difffcrcntial geometry, since in three-dimensions, 
the shadow space K is one-dimensional, so the sectional curvature (regarded as an 
cndomorphism of K) is automatically proportional to the identity endomorphism, 
so is pure trace. So non-trivial Weyl curvature arises first in dimension four, just 
as in metric differential geometry. 

B.5. Example: indefinite metrics. We consider the case that G — 2~ 1 g ab (x)v a v b , 
where g a b{x) = gbaix) and the metric g a b{x) is invertible and indefinite; these com- 
prise the Kleinian (ultra-hyperbolic) metrics and the Lorentzian (hyperbolic) case 
of general relativity. Then: 

g a = g a bv b , g a b = g a b, g a bc = 0, G a = 2~ 1 v b v c d a gbc, 
u a g a b = G b - v c d c g b = 2~ 1 v a v c d b g ac - v a v c d c g ab 
= 2- 1 v c v d {d b9cd - 2d (c9d)b ) = -v c v d T cd a g ab , 
Y bc a = -2-y d (9 d5fcc - 2d (b g c)d ). 
So T cd a are the usual Christoffcl symbols of the Levi-Civita connection. Then: 

u a = -v b v c r bc a , u b a = -±D b u a - v c r bc a , 

V = v a d a - v b v c T bc a D a , H a = a - v c T ab c D c . 
So V is the ordinary null geodesic spray, whereas the vector fields H a represent the 
horizontal vector fields on the tangent bundle of the standard Levi-Civita connec- 
tion. A conventional conformal transformation has J(x,v) = J(x) ^ 0, so j a = 
and j a b = 0. Also p = k = 2 and q = and we have: 

H = -h ab v a v b = GJ(x), h ab = J(x)g ab , 

u a = -f a + G( a + G 2 C + G 3 <r a + (d + Gr] + G 2 P )v a 1 
7 a = h ab (H b - v c d c h b ), W = v a 8 a + l a D a , 5 = J- x v a J a , 

C = -h ab j b , r = ^ a = o, v = p = o, 

u a = h ab (H b - v c d c h b ) - HJ- l h ab J b + v a J- 1 v b J b . 
Write A bc a for the Christoffel symbols of h ab . Then we have: 

-T bc a v b v c = v»v c (-A bc a - 2- 1 h bc J- 1 h ad J d + SfrJ- 1 J c) ), 

Y bc a = A bc a + -J^ 1 h bc h ad J d - J _1 J((,5^. 
This agrees with the standard transformation law for Christoffel symbols: 

A bc a = -2- x h ad {d A h bc - 2d (b h c)d ) = -2~ 1 J~ 1 g ad (d d (Jg bc ) - 2d (b (Jg c)d )) 
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= r 6c a -^J 1 g ad (gb C Jd-2J(b9c)d) = r fcc a - 1 h ad h bc J d + J 1 J( b 5* ) . 

B.6. The Lie derivatives of A. We have in local co-ordinates: 

A = dx a ® D a . 

Then we have: 

££ v \ = dv a ®D a + dx a ® (-d a - D a u b )D b , 
I = dv a ® D a + dx a <g> d a , 

r = 1(J + jSf v A) = (dv a - ^dx b D b u a ) ®D a = ^ a ®D a , 

cj) a = dv a + Ujfdx b , 

t£vY= X -^\ = % v (4, a ®D a ) 

= L V (dct> a ) ®D a -cj> a ® (8 a - 2U h a D b ) 
= V(U£)dx b ®D a - v b (dU£) ®D a -4> a ® (da - 2U b D b ) 
= V(U£)dx b ®D a - d(v b U£) ®D a + U£(dv b ) ®D a -4> a ® (H a - U b D b ) 
= V(U b a )dx b ®D a + (du a ) ®D a + U^(4> b - U b dx c ) ® D a - (j> a <8> (H a - U b D b ) 
= V(U b a )dx b ®D a -2^ b (U b 1 )®D a +dx b (H b u a )®D a +UZ(<t> b ^ 

= (V(Ut) + H b u a - U?U§)dx b ®D a -4> a ®H a 
= T£dx b ®D a ~cf> a <E>H a . 

So we have: 

j2Vr = = T£dx b ® £>„ - 0° ® H a . 
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